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LIST  OF  VARIABLES  USED  IN  P-STAT 


NSECT: 

LA: 

LANT: 

NCOUP: 

ESTO: 

WSTO: 

RSTO: 

NRUN: 

lOFF: 

IT: 

XN: 

SPD: 

ALT: 

MODEF: 

FSTRT: 

FSTP: 

FDEL: 


An  integer  variable  specifying  the  program  option  to  be 
used  (corona  noise  or  streamer  noise) . 

An  integer  variable  specifying  the  antenna  location. 

An  alphanumeric  variable  describing  the  antenna  location. 

An  integer  variable  specifying  the  number  of  coupling 
coefficients  to  be  read  from  data  cards. 

A  floating  point  array  containing  the  NCOUP  antenna- 
elevator  coupling  coefficients. 

A  floating  point  array  containing  the  NCOUP  antenna-wing 
coupling  coefficients. 

A  floating  point  array  containing  the  NCOUP  antenna-rudder 
coupling  coefficients. 

An  integer  variable  specifying  the  number  of  program 
cycles  to  be  made  using  the  same  coupling  coefficients. 

An  integer  variable  specifying  the  locations  of  the  p -static 
discharges  which  are  to  be  considered  "quiet". 

An  alphanumeric  variable  describing  the  type  of  aircraft 
under  investigation. 

A  floating  point  variable  specifying  the  size  of  the 
aircraft  relative  to  a  KC-135. 

A  floating  point  variable  specifying  the  aircraft  speed. 

A  floating  point  variable  specifying  the  aircraft  altitude. 

An  integer  variable  specifying  the  frequency  select  mode 
the  user  wishes  to  use  (uniform  or  non-uniform  frequency 
Intervals) . 

(If  MODEF  equals  0)  A  floating  point  variable  specifying 
the  desired  starting  frequency  (in  MHz). 

(If  MODEF  equals  0)  A  floating  point  variable  specifying 
the  desired  stopping  frequency  (in  MHz) . 

(If  MODEF  equals  0)  A  floating  point  variable  specifying 
the  frequency  increment  between  FSTRT  and  FSTP  (in  MHz). 
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NFR;  (If  MODEF  does  not  equal  0)  An  integer  variable  specifying 

the  number  of  user-selected  frequencies  to  be  read  in  from 
cards . 

FREQU:  (If  MODEF  does  not  equal  0)  A  floating  point  variable 

specifying  the  user-selected  frequency  (in  MHz).  The 
maximiua  number  of  F^QU  cards  allowed  is  90. 

AANT:  A  floating  point  variable  specifying  the  antenna  induction 

area  (in  square  meters). 

BNDW:  A  floating  point  variable  specifying  the  receiver  bandwidth 

(in  kHz) . 

ICLO;  An  integer  variable  specifying  the  type  of  cloud  the 
aircraft  is  flying  through. 

IC:  An  alphanumeric  variable  describing  the  type  of  cloud  the 

aircraft  is  flying  through. 

(Variables  Used  Only  in  Streamer-Noise  Calculations) 

IM:  An  integer  variable  specifying  the  type  of  dielectric 

material  being  charged. 

IMAT:  An  alphanumeric  variable  describing  the  type  of  dielectric 

material  being  charged, 

DAFT:  A  floating  point  variable  specifying  the  distance  (in  meters) 

the  receiving  antenna  is  located  behind  the  windshield 
canopy  or  the  radome. 

WX:  A  floating  point  variable  specifying  the  minimum  characteristic 

dimension  (in  meters)  of  the  dielectric  material  being  charged. 

DIERAT;  A  floating  point  variable  specifying  the  ratio  of  the  frontal 
area  of  the  dielectric  material  to  the  frontal  area  of  the 
aircraft . 
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I  INTRODUCTION 


When  an  aircraft  or  other  flight  vehicle  is  operated  in  precipita¬ 
tion  containing  ice  crystals  or  other  particulate  materials,  frictional 
electrification  associated  with  particle  impact  causes  the  impinging 
particles  to  acquire  a  net  charge  and  to  deposit  an  equal  and  opposite 
charge  on  the  vehicle.^" The  charging  occurs  on  the  frontal  metallic 
and  dielectric  portions  of  the  vehicle.®’*^  Although  the  charge  deposited 
by  a  single  ice  crystal  changes  the  potential  of  the  aircraft  only 
slightly  (of  the  order  of  0.01  volt  for  the  case  of  a  KC-135  struck  by 
a  cirrus-cloud  crystal),^  the  particle  impact  rate  in  a  typical  cloud 
is  sufficient  to  cause  the  vehicle  potential  to  reach  hundreds  of  kilo¬ 
volts  in  less  than  a  second.^ 

The  electrification  of  the  vehicle  is  of  relatively  little  concern 
in  itself  because  the  energies  involved  are  small,  and  since  the  electro¬ 
static  fields  do  not  penetrate  to  the  interior.  It  is  the  consequences 
of  the  electrification  that  are  of  concern  to  the  EMC  engineer.  When 
the  vehicle  potential  reaches  roughly  100  kV,  the  electric-field  intensity 
at  the  aircraft  extremities  becomes  sufficiently  high  that  electrical 
breakdown  of  the  air  (corona  discharge)  occurs.®  At  aircraft  operating 
altitudes,  the  corona  breakdown  from  the  extremities  occurs  not  as  a 
continuous  flow  of  charge,  but  as  a  series  of  pulses  with  roughly  10  ns 
rise  times  and  200  ns  duration  and  therefore  generates  radio  noise  over 
a  broad  spectrum.^ ’ ® 


References  are  listed  at  the  end  of  this  Users  Manual. 
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A  similar  situation  exists  on  the  dielectric  frontal  surfaces. 

As  charge  continues  to  accumulate  on  the  dielectric,  the  potential  to 
the  airframe  rises  until  the  electric-field  intensity  at  the  dielectric 
surface  becomes  sufficiently  high  that  voltage  breakdown  (streamer 
discharge)  across  the  plastic  surface  occurs.  A  surface  streamer 
involves  the  rapid  transfer  of  charge  over  a  substantial  distance, 
and  also  generates  serious  radio  frequency  interference.®’ 

The  degree  to  vdiich  the  radio  frequency  noise  generated  by  corona 
and  streamer  discharges  couples  into  electronic  systems  on  the  flight 
vehicle  is  determined  by  the  relative  locations  of  the  noise  source  and 
the  ''antenna”  via  which  the  noise  is  coupled  into  the  affected  system. 

In  addition,  the  coupling  depends  upon  frequency,  the  size  of  the 
vehicle,  and  the  size  of  the  antentia."^’ 

On  earlier  efforts,  various  aspects  of  the  problem  of  precipitation- 
static  noise  generation  and  coupling  were  studied  analytically  and 
experimentally  both  in  the  laboratory  and  in  flight  tests.  Unfortunately, 
the  results  of  these  efforts  are  spread  over  a  large  number  of  reports, 
each  of  which  treats  only  a  limited  part  of  the  overall  problem.  Thus 
the  EMC  engineer  is  in  the  position  of  having  to  be  familiar  with  all 
of  the  publications  in  considerable  depth  if  he  is  to  apply  the  results 
of  the  earlier  work  to  his  problems. 

In  order  to  overcome  these  problems,  SRI  developed  a  computer 
program,  entitled  PSTAT,  which  will  accurately  predict  the  effects  of 
p-static  noise  in  aircraft  systems.  The  computer  program  has  been 
demonstrated  to  allow  the  EMC  engineer,  or  systems  designer,  to  determine 
the  effects  of  p-static  charging  on  a  wide  variety  of  aircraft  types  and 
under  a  wide  variety  of  flight  regimes.  Since  the  program  is  based  on 
the  results  of  earlier  experimental  and  analytical  work,  the  limitations 
of  this  earlier  work  are  reflected  in  the  computer  program.  The  accuracy 
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of  PSTAT  depends  on  the  modelling  and  on  the  faithfulness  with  which  the 
experimental  analytical  data  represent  the  true  picture  of  p-static  noise* 
It  is  felt  that  PSTAT  is  accurate  to  within  a  few  percent  for  KC-135 
type  aircraft,  decreasing  to  tens  of  percent  for  widely  divergent  air¬ 
craft  types  (delta  wing  fighters,  for  example).  Although  it  has  been 
possible  to  extend  the  applicability  of  the  first-generation  program 
described  here  somewhat  beyond  the  strict  confines  of  the  earlier  work, 
there  are  situations  in  which  the  program  simply  cannot  be  applied.  For 
example,  with  the  present  program,  it  is  not  possible  to  consider  heli¬ 
copters  or  rockets  because  their  geometries  are  radically  different  from 
aircraft. 

This  users  manual  is  intended  to  guide  the  program  user  through  the 
input  and  output  requirements  of  the  program.  Sample  input  decks  and 
output  listings  are  included  in  this  users  manual  to  help  the  user 
understand  the  proper  input-deck  setup.  Specific  modeling  techniques 
are  not  explained  in  this  manual  because  they  are  fully  explained  in 
the  accompanying  Final  Report  under  this  contract. 

The  philosophy  applied  in  creating  the  present  program  was  one  of 
simplicity.  The  authors  felt  that  direct  in-line  coding  was  more  ap¬ 
propriate  to  the  needs  of  potential  users  than  were  more  complicated 
coding  techniques.  In-line  coding  affords  the  non-programmer  user  the 
convenience  of  being  able  to  look  at  the  program  and  determine  the 
sequence  of  events  that  have  just  taken  place  and  those  that  are  about 
to  begin. 

Extensive  comments  have  been  inserted  throughout  the  program  in 
order  to  clarify  the  various  program  steps. 
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II  HARDWARE  REQUIREMENTS  AND  LANGUAGE 

A*  Hardware  Requirements 

PSTAT  was  designed  to  run  with  a  minimum  computer  configuration* 
The  program  uses  a  card  reader  for  input  and  a  line  printer  for  output. 
No  additional  peripherals  are  required. 

The  program  uses  5203^^^  words  of  core  storage. 

Execution  time  is  dependent  on  the  parameters  selected  during 
input,  but  typical  execution  times  of,  perhaps,  5  to  10  seconds  could 
be  expected  for  typical  calculations,  and  this  time  would  include  the 
card  read,  CPU,  and  printer  times. 

It  is  estimated  that  the  CPU  time  required  for  a  typical  run  is 
on  the  order  of  100  ms. 

B.  Language 

PSTAT  is  written  in  standard  ASA  FORTRAN. 
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Ill  COMPUTER  PROGRAM 


A.  General 

The  experimental  and  analytical  data  regarding  p-static  noise  is 
discussed  fully  in  Section  II  of  the  Final  Report  (Vol.  1)  written  under 
this  contract  and  will  not  be  repeated  here. 

The  nature  of  the  material  presented  in  the  final  report  was  such 
that,  in  some  cases,  exact  analytical  expressions  could  be  used  in  the 
computer  program.  In  other  cases,  approximations  to  the  desired  param¬ 
eters  were  used;  and  in  still  others,  where  the  data  did  not  lend  them¬ 
selves  to  approximation,  the  data  were  simply  stored  in  tabular  form. 

B.  Flowchart 

Based  on  immediate  needs,  the  requirements  anticipated  in  the  future, 
and  information  currently  available,  a  flowchart  was  developed  to  be  a 
guideline  for  the  coding  effort.  This  flowchart  is  shown  in  Figure  1. 

It  can  be  seen  from  this  figure  that  the  p-static  program  is  broken 
into  two  sections,  or  modules.  Module  1  deals  with  the  calculation  of 
noise  generated  in  antennas  by  corona  discharges  from  the  aircraft 
extremities.  Module  2  deals  with  the  calculation  of  noise  generated  in 
antennas  by  surface  streamer  discharges  across  the  plastic  surfaces  of 
the  aircraft's  radomes  and  canopies. 

During  program  execution,  either  Module  1  or  Module  2  is  selected 
by  the  user  by  use  of  a  data  card  read  in  as  the  first  data  card. 

It  can  be  observed  from  this  figure  that  an  input  data  error  test 
is  made  only  on  the  data  input  to  Module  1.  It  was  decided  that  the 
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FIGURE  1  PSTAT  FLOWCHART 
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input  requirements  of  Module  2  were  sufficiently  simple  that  an  input 
error  check  could  not  be  justified,  whereas  the  input  requirements  of 
Module  1,  while  not  complex,  were  sufficiently  confusing  to  warrant 
the  error  check* 

A  brief  description  of  the  contents  of  each  program  module  is 
given  below.  The  input  and  output  details  of  each  module  are  not 
discussed  here,  but  are  left  for  a  later  section  of  this  manual.  The 
mathematical  processes  of  the  calculations  performed  in  the  modules  are 
fully  described  in  the  Final  Report,  so  they  will  not  be  repeated  here. 


C.  Module  1 — Corona  Noise 
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After  the  data  cards  have  been  input,  an  error  check  is  made  on 
several  of  the  important  parameters  of  the  program.  PSTAT  will  produce 
the  error  message 


INPUT  ERROR**** 

print  the  input  deck,  repeat  the  error  message,  and  then  halt,  if  any 
of  the  following  errors  are  detected: 

•  More  than  100  coupling  coefficients  for  each  extremity 
are  either  read  into  the  program  or  requested  to  be 
read  into  the  program. 

•  More  than  90  frequencies  have  been  read  into  the  program 
or  requested  to  be  read  into  the  program  (for  MODEF  .NE.O). 
(Note:  For  MODEF  .EQ.O  any  number  of  frequencies  may  be 
evaluated — see  description  of  constants  and  variables 
below.) 

•  The  requested  frequency  ranges  and/or  frequency  interval 
are  not  consistent--e.g. ,  if  the  last  frequency  were 
smaller  than  the  first  frequency,  or  if  Af  were  0  or 
negative — note:  (This  check  is  made  only  if  MDDEF  .EQ.O). 

•  The  discharge  quench  code  does  not  reflect  any  of  the 
possible  quenching  modes. 

•  The  aircraft's  altitude  is  greater  than  80,000  ft. 
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After  the  input  deck  has  passed  the  error  check,  it  is  printed 
out,  showing  the  user  the  parameters  he  has  selected  for  evaluation. 

The  next  step  in  the  program  is  the  calculation  of  the  total 
charging  current  to  the  aircraft.  The  total  charging  current  (and 
hence  discharging  current  in  the  steady-state  case  under  consideration) 
is  determined  from  the  aircraft  speed,  its  size  (relative  to  a  KC-135), 
and  the  type  of  cloud  it  is  penetrating.  At  the  same  time  the  charging 
current  is  calculated,  the  probability  that  this  charging  current  will 
be  exceeded  is  also  calculated. 

Since  the  noise  coupled  into  the  antenna  is  a  function  of  the 
antenna  induction  area  and  aircraft  size,  the  coupling  coefficients 
are  then  scaled  to  reflect  the  antenna  induction  area  and  the  aircraft 
size.  The  next  step  in  the  program  distributes  the  total  charging 
current  among  the  extremities  (rudder,  elevator  tips,  wing-tips)  and 
then  calculates  the  discharge  source  spectrum  normalizers,  \diich  are 
used  to  determine  the  intensity  of  the  corona  spectra. 

After  the  pressure  (altitude)  and  frequency  parameters  have  been 
initialized,  the  equivalent  noise-field  calculations  begin.  The  spectral 
function,  PREL,  in  the  program  is  calculated  using  the  approximations 
detailed  in  the  Final  Report,  and  the  coupling  data  are  linearly 
interpolated  from  the  table  of  coupling  coefficients  established  during 
the  input  phase  of  the  program. 

After  the  short-circuit  antenna  current  and  equivalent  noise  fields 
have  been  calculated  they  are  printed  out  for  the  frequency  currently 
being  investigated.  A  frequency- increment  test  directs  the  program 
either  to  a  "continue  processing"  statement  or  to  a  "completion" 
statement. 
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D,  Module  2 — Streamer  Noise 

The  technique  used  to  calculate  the  equivalent  noise  caused  by 
streamer  discharge  closely  parallels  the  technique  used  to  calculate 
corona  noise.  After  the  data  cards  have  been  read  in,  the  input  deck  is 
printed  showing  the  user  the  parameters  he  has  selected  for  evaluation. 
This  serves  as  an  error  check  on  the  input  data. 

The  next  step  in  the  program  is  the  calculation  of  the  total 
charging  current  to  the  aircraft.  The  total  charging  current  is  deter¬ 
mined  from  the  aircraft  speed,  size,  and  type  of  cloud  it  is  penetrating. 
At  this  same  time,  the  probability  that  this  charging  current  will  be 
exceeded  is  also  calculated. 

The  next  step  in  the  program  is  the  calculation  of  the  streamering 
current.  The  streamering  current  is  given  by  the  ratio  of  the  dielectric 
surface  frontal  area  to  the  total  aircraft  frontal  area  multiplied  by 
the  calculated  aircraft  charging  current. 

After  the  frequency  parameters  have  been  initialized,  the  streamer 
spectrum  is  calculated  at  the  particular  frequency  being  examined.  The 
short  circuit  antenna  current  is  then  calculated  and  the  equivalent 
noise  field  is  finally  obtained  and  printed  out.  A  frequency  increment- 
test  directs  the  program  either  to  a  "continue"  processing,  or  a 
"completion"  statement. 

The  inherent  qualities  of  program  PSTAT  are  that,  in  the  brief 
module  descriptions  given  above,  many  years  of  accumulated  experimental 
data  have  been  combined  to  form  a  unified  program  to  solve  many  types 
of  problems  involving  precipitation-static-induced  noise  in  avionics 
systems.  While  the  program,  taken  in  its  entirety  involves  considerable 
sophistication,  the  individual  calculations  are  quite  simple  and  easily 
followed  in  the  program  documentation.  Accordingly,  we  have  not  provided 
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flow  charts  for  the  calculation  of  every  parameter  because  it  was  felt 
that  they  would  be  simple  but  so  numerous  as  to  detract  from  the  utility 
of  this  manual. 
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IV  INPUT 


PSTAT  utilizes  three  input  areas:  (1)  The  initial  one-card  input 
to  specify  Module  1  (corona  noise)  or  Module  2  (streamer  noise),  (2)  the 
input  area  for  the  corona-noise  calculation,  and  (3)  the  input  area  for 
the  streamer-noise  calculation. 

At  any  one  time  the  user  will  use  only  two  of  these  areas:  The 
module-select  area  and  the  corona-noise  area,  or  the  module- select  area 
and  the  streamer-noise  area. 

The  requirements  and  formats  for  each  of  these  areas  are  given 
below.  The  order  in  which  the  material  is  presented  is  the  order  in 
\diich  the  input  deck  should  be  arranged. 

A.  Module  Select  Area 

•  Card  1 — This  will  always  be  the  first  card  of  the  data  deck, 
and  it  contains  either  a  1  (Module  1),  or  a  2  (Module  2)  and 
directs  the  program  to  the  desired  module.  The  card  should 
be  in  an  II  format. 


B.  Corona-Noise  Module 

The  description  of  each  of  the  cards  to  be  input  into  this  module 
is  given  below,  in  the  order  of  their  location  in  the  input  deck. 

•  Card  2— LA,  LANT;  Format  II,  IX,  7A2 

LANT  is  a  14-character  alphanumeric  briefly  describing  the 
location  of  the  antenna  under  test  (i.e.,  BELLY,  FUSELAGE, 
TAILCAP,  etc.)  and  is  used  only  for  output  annotation. 

LA  is  a  single-digit  fixed-point  variable  describing  the 
antenna  location.  Set  LA  =  0  if  the  antenna  is  not  located 
at,  or  near,  an  extremity  (e.g.,  a  belly -mounted  antenna). 
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If  the  antenna  is  located  at,  or  near,  the  elevator 
extremity,  set  LA  =  !•  If  the  antenna  is  located  at, 
or  near,  a  wing-tip,  set  LA  =  2.  Set  LA  =  3  if  the 
antenna  is  located  at,  or  near,  the  rudder  extremity* 

This  parameter  is  used  to  scale  the  coupling  coef¬ 
ficients  to  the  scale  size  of  the  aircraft,  for  those 
discharge  locations  not  located  near  the  antenna*  The 
coupling  coefficient  describing  the  coupling  between 
noise  sources  and  extremity-located  antennas  is  not 
scaled  to  aircraft  scale  size  if  the  antenna  is  located 
near  those  noise  sources*  The  other  coupling  coefficients, 
however,  are  scaled,  and  the  reasons  for  scaling  are 
described  in  the  final  report. 

•  Card  3— NCOUP;  Format  13 

This  is  a  fixed-point  number  specifying  the  number  of 
coupling  coefficients  to  be  read  from  cards  (Maximum  =  100)* 

•  Card  4--EST0,  WSTO,  RSTO;  Format  3(E9.2,1X) 

These  are  the  array  names  for  the  storage  of  the  NCOUP 
coupling  coefficients.  The  data  on  these  cards  are 
experimentally  derived  quantities  and  until  the  user  gains 
familiarity  with  the  program,  or  until  more  data  become 
available,  the  SRI-supplied  decks  of  coupling  coefficients 
should  be  used.  The  user  should  note  that  SRI  has  supplied 
two  decks  of  coupling  coefficients:  one  for  extremity-to- 
tail-cap  antennas;  and  one  for  extremity-to-belly  antennas* 
The  user  should  select  the  deck  appropriate  to  his  needs — 
tail-cap  or  belly-mounted  (fuselage-mounted)  antennas* 

•  Card  5--NRUN;  Format  13 

This  card  specifies  the  number  of  program  cycles  to  be  made 
using  the  same  coupling  data  but  various  other  parameters. 

It  is  suggested  that  until  the  user  is  familiar  with  the 
program,  NRUN  be  limited  to  1. 

•  Card  6--I0FF;  Format  II 

This  card  specifies  which  (if  any)  of  the  corona  discharges 
should  be  suppressed  by  40  dB.  (40  dB  is  typical  of  the 
quieting  provided  by  p-static  dischargers  on  aircraft.) 

The  codes  are  as  follows: 
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lOFF  =  1  All  discharges  permitted 

lOFF  =  2  Rudder  discharge  quieted  by  40  dB 

lOFF  =  3  Wing-tip  discharges  quieted  by  40  dB 

lOFF  =  4  Elevator-tip  discharges  quieted  by  40  dB 

lOFF  =  5  Rudder  and  wing-tip  discharges  quieted  by  40  dB 

lOFF  =  6  Rudder  and  elevator-tip  discharges  quieted  by 

40  dB 

lOFF  =  7  Elevator  and  wing-tip  discharges  quieted  by 
40  dB. 

•  Card  7--IT;  Format  6A2 

This  is  a  12-character  alphanumeric  describing  the  type  of 
aircraft  under  investigation  (i.e.,  TRANSPORT,  FIGHTER,  etc.), 
and  is  used  only  for  output  annotation. 

•  Card  8--XN,  SPD,  ALT;  Format  F5.2,  IX,  F6.1,  IX,  F4.1 

This  card  contains  the  information  about  the  aircraft’s 
size,  XN  (relative  to  a  KC-135),  and  its  speed  (in  mph) 
and  its  operating  altitude  (in  kft). 

•  Card  9— MODEF;  Format  II 

This  card  specifies  the  frequency-select  mode  the  user  wishes 
to  use.  If  MODEF  equals  0,  it  means  that  the  user  has  decided 
to  use  uniformly  spaced  frequency  intervals.  If  MODEF  is  not 
equal  to  0,  it  means  that  the  user  has  decided  to  use  frequencies 
that  will  be  read  in  from  cards  at  a  nonuniform  Af. 

•  Card  10— (If  MODEF  .EQ.O)  FSTRT,  FSTP,  FDEL;  Format  3(F5.2,  IX) 

This  card  contains  the  desired  starting  frequency  (in  MHz), 
ending  frequency  (in  MHz),  and  frequency  increment  (in  MHz) 
if  MODEF  is  equal  to  zero. 

•  Card  10— (If  MODEF  .NE.O)  NFR;  Format  13 

This  card  specifies  the  number  of  user-selected  frequencies  to 
be  read  into  the  program.  (The  maximum  number  allowed  is  90.) 

•  Cards  10a,  10b,  10c,  etc. --(If  MODEF  .NE.O)  FREQU;  Format  E9.2 

These  cards  are  the  user-selected  frequencies  (in  MHz). 

There  should  be  NFR  of  these  cards. 
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•  Card  11--AANT,  BNDW;  Format  2(F5.2,  2X) 

This  card  contains  the  information  specifying,  the  receiving 
antenna's  induction  area  (in  m^)  and  the  receiver  bandwidth 
(in  kHz) . 

•  Card  12— ICLO,  IC;  Format  II,  IX,  7A2 

This  card  contains  the  information  about  the  type  of  particulate 
material  the  aircraft  is  flying  in. 

ICLO  =  1  implies  a  cirrus  cloud  or  low  charging  material. 

ICLO  =  2  implies  a  stratocumulus  cloud  or  moderate 
charging  material. 

ICLO  =  4  implies  a  snow  cloud  or  high-charging  material. 

IC  is  a  14-character  alphanumeric  description  of  the  cloud 
material.  It  is  used  only  for  output  annotation. 

C.  Streamer-Noise  Module 

•  Card  2--LANT;  Format  4A2 

This  alphanumeric  is  described  in  Section  IV-B  above. 

•  Card  3- -IT;  Format  6A2 

This  alphanumeric  is  described  in  Section  IV-B  above. 

•  Card  4--XN,  SPD,  ALT;  Format  F5.2,  IX,  F6.1,  IX,  F4.1 

The  data  on  this  card  are  described  in  Section  IV-B  above. 

•  Card  5--M0DEF;  Format  II 

The  data  on  this  card  are  described  in  Section  IV-B  above. 

•  Card  6-- (If  MODEF  .EQ.O)  FSTRT,  FSTP,  FDEL;  Format  3(F5.2,  IX) 
The  data  on  this  card  are  described  in  Section  IV-B  above. 

•  Card  6— (If  MODEF  .NE.O)  NFR;  Format  13 

The  data  on  this  card  are  described  in  Section  IV-B  above. 

•  Card  6a,  6b,  6c— (If  MODEF  .NE.O)  FREQU;  Format  E9.2 

The  data  on  these  cards  are  described  in  Section  IV-B  above. 
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•  Card?— AMT,  BNDW;  Format  2(F5.2,  2X) 


The  data  on  this  card  are  described  in  Section  IV-B  above. 

•  Card_8.-ICL0,  IC;  Format  II,  IX,  7A2 

The  data  on  this  card  are  described  in  Section  IV-B  above. 

•  Card  9— IM,  IMAT;  Format  II,  IX,  7A2 

This  card  contains  the  information  about  the  type  of 
dielectric  material  being  charged. 

IM  =  1  implies  that  a  windshield  (canopy)  is  being 
charged. 

IM  =  2  implies  that  a  radome  is  being  charged. 

IMAT  is  a  14-character  alphanumeric  description  of  the 
dielectric  material  (i.e.,  WINDSHIELD,  or  RADOME).  It  is 
used  only  for  output  annotation. 

•  Card  10— DAFT, WX;  Format  2 (F5. 2,  2X) 

This  card  describes  the  antenna  location  with  respect  to 
the  charging  material,  and  the  minimum  characteristic  dimen¬ 
sion  of  the  dielectric  material  being  charged. 

DAFT  specifies  the  distance  (in  meters)  the  receiving  antenna 
is  located  behind  the  windshield  canopy  or  the  radome.  If 
the  receiving  antenna  is  located  immediately  beneath  the 
dielectric  material,  DAFT  should  be  read  in  as  0.00  m. 

WX  specifies  the  minimum  characteristic  dimension  (in  meters) 
of  the  dielectric  material  being  charged- -i.e.,  the  width  of 
a  rectangular  section  of  dielectric.  The  floating-point 
variable,  WX,  may  be  thought  of  as  roughly  twice  the  length 
of  the  longest  possible  streamer  discharge  on  the  dielectric 
region  under  consideration. 

•  Card  11--DIERAT;  Format  F5.2 

DIERAT  is  the  ratio  of  the  frontal  area  of  the  dielectric  to 
the  frontal  area  of  the  aircraft. 

In  the  event  that  windshield  canopy  streamering  is  being 
considered,  DIERAT  should  specify  the  ratio  of  the  total 
frontal  area  of  the  dielectric  to  the  total  frontal  area 
of  the  aircraft. 
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If  radome  streamering  is  being  considered,  DIERAT  should 
specify  the  ratio  of  the  radome *s  forward  3  feet  of  area 
to  the  total  frontal  area  of  the  aircraft. 

It  can  be  seen  from  the  input  requirements  described  above  that 
the  use  of  alphanumerics  has  been  limited  to  annotation  only,  while 
parameters  which  affect  the  processing  has  been  limited  to  BCD  (numbers). 
This  technique  could  have  been  changed  so  that  alphanumerics  directed 
some  of  the  processing,  but  it  was  felt  that  this  would  confuse  the 
input  requirements  of  PSTAT.  The  example  INPUT/OUTPUT  shown  later  in 
this  volume  will  illustrate  the  use  of  the  BCD/Alphanumerics  input 
data  described  above. 
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V  OUTPUT 

During  output,  the  user-supplied  quantities  that  affect  the  computed 
results  are  printed  out  before  the  induced  equivalent  noise  fields  are 
printed  out. 

If  an  error  is  detected  during  the  processing  of  the  corona-noise 
input  deck,  an  error  message  is  produced.  No  error  checks  are  made 
during  the  processing  of  the  streamer-noise  input  deck,  since  the  input 
requirements  for  this  module  are  quite  simple. 

After  the  input  quantities  have  been  listed,  the  charging  current 
is  calculated  and  printed  out.  The  probability  that  the  charging  current 
will  exceed  the  calculated  value  (for  the  specified  conditions  of  altitude, 
speed,  aircraft  size,  and  cloud  type)  is  also  calculated  and  printed  out. 

The  short-circuit  currents  induced  in  the  receiving  antenna  and  the 
associated  equivalent  noise  fields  are  then  calculated  and  printed  out 
for  all  of  the  user-desired  frequencies.  The  dimensions  of  these  output 
quantities  are  megahertz  and  hertz  for  the  user-specified  frequencies, 
amperes  for  the  short-circuit  current,  and  volts  per  meter  for  the 
equivalent  noise  fields. 

It  should  be  noted  here  that  if  the  user  elects  to  use  the  streamer- 
ing  model  for  an  antenna  immediately  beneath  the  canopy  or  radome,  no 
equivalent  noise  field  is  calculated  or  printed.  The  reasons  for  this 
are  fully  described  in  the  final  report. 

Examples  of  the  output  are  given  in  a  later  section  of  this  manual. 
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VI  SAMPLE  INFUT/OUTPUT 


This  section  gives  several  examples  of  the  use  of  program  PSTAT, 
together  with  example  input  deck  setup  and  output  listing. 

A«  Example  1 

Calculate  the  equivalent  noise  field  induced  in  an  antenna  on  the 

tail-cap  of  a  KC-135  transport  aircraft.  Assume  that  the  antenna  has 

2 

an  induction  area  of  8.6  m  ,  and  that  the  receiver  has  a  bandwidth  of 
1  kHz.  -  Further  assume  that  the  aircraft  is  flying  at  a  speed  of  600  mph 
at  an  altitude  of  20,000  feet  through  cirrus  cloud.  Allow  all  extremities 
of  the  aircraft  to  discharge  and  evaluate  the  equivalent  noise  fields 
at  uniformly  spaced  frequencies  of  from  0.1  MHz  to  4.0  MHz  in  steps 
of  0.1  MHz. 
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1.  Input  Deck 


The  input  deck  required  to  evaluate  this  problem  is  as 


follows : 


1 


3«TAILCAP 

15* 


+0.A1E-03 

+0.?3E-03 

+  0 

♦35E-C1 

0 

TAILCAP 

+0.35E-03 

+0.3CE-0? 

+  0 

• "C 1 

1  ^'HZ 

TAICCAP 

+C.2CE-C3 

+0 • 56E«03 

+  C 

.35E-C1 

a 

TAILCAP 

♦0.3CE-02 

+0.16E-0? 

.SSE-Cl 

2  ►"•HZ 

tailcaf 

.  +0.50E-02 

*  35t  *01 

4  MHZ 

TAILCAP 

+0.27E-03 

+0.11E-02 

+c 

.372-01 

E  '^MZ 

tailcaf 

♦0.27E-D2 

+0#75E-C2 

+ 1 

.403-01 

6 

’’ailcap 

+C.32E-02 

+  C.tlCE*C? 

+  r 

.392-01 

7  HW2 

tailcaf 

+0.A3E-02 

+0.i7E-0? 

4-2 

.??F.-ci 

f  MHZ 

tailcaf 

+C.70E-C2 

-»'0«icE-oa 

+  r 

.353-01 

9  ►"HZ 

tailoaf 

+0.1CE-01 

♦0. ACE”"3 

4-  C 

. 333-01 

IC^'HZ 

tAILCAP 

♦O.ISE-Ol 

4  D 

.402-01 

ll'"HZ 

^AILCAP 

+0* 13E"C1 

+6.742-03 

4  r 

.  f  1 2  -  0 1 

12^^HZ 

■"AIlCAP 

+0.12E-01 

+0.3CE-C-3 

+: 

.571-01 

13^"wz 

TAILCAP 

+0»lbE-bi 

+0. i CE-0? 

4  0 

.553-02 

14MHZ 

TAILCAP 

h _ _ 

1 

KC-135 

1.00  6cb.'c  rC.D 

o’.io  4.0C  c.ic 
8.6  *  1.0 
1-CIRPUS  CL?U0 
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BEST  available  COPY 


2.  Output  Deck 


The  program  output  is  as  follows: 


S?±.l3ltllQ^£Lt£JRl£jn..mSi.lL 


P-STATIC  E'/Al^UATED  PBR  A  KC-13?  AIRCRAFT 

WITH  THE  RECEIVINj  ANTEMHA  ueCATEC  AT  THE  TAILCAP 


SCALE  SIZE 

SPEED 

ALTITUDE 

CLOUD  TYPE 

ri^PHD 

C‘<FT3 

1*00 

600*0 

20*0 

CIRRUS  cloud 

start  FRECi 

STOP  EREG* . 

DELTA-F 

. [HHZ  ]  , 

r^'HZ] 

.10 

4 .  DC 

.10 

RECEIVER 

ANTENNA 

HQfSE  “  . . 

INDUCTION 

BANDWIDTH 

APE  A 

tKHZl 

CM**2I 

. . 

. 8.6C 

ALL  discharges 

PE^^ITTED 

THr  CALC'JLATrr  C-'At;3INS  CURR^M  TS  l«000E-03  AMPS  .. 

THE  »R9B_A3ILITX.  IS  __jlOCSC  THjy:  THE  CHARGING  CURPENT 
WILL  3E  SREATE;?  THAN  i.CGCF-Oa  AHRS 


FRES'JENC'^ 

■SHf^P 

T-CIRCyiT 

EOyiVALENT 

EQUIVALENT 

CU^^ENT 

NeiSE  FIELD 

N6ISE  field 

CNHZI 

.IHZ3 

c 

A.N®S,3 

tV8LTS/M.] 

CDBV/M3 

•  1C 

1 •CCCE 

r.P 

g 

.434E-C7 

1.765E-05 

-3.506E  01 

.?o 

5.0CCE 

05 

g 

.376E-C7 

S.765E-03 

-4.114E  01 

•  30 

3.0CCE 

GE 

g 

.5P1E-C7 

5.777E-03 

-4.476E  01 

•  40 

4. 'COE 

05 

fi 

.154E-C7 

4.267E-03 

-4,739E  01 

.50 

.  S*COCE 

w  -■ 

7 

.9S9E-C7 

3.348E-03 

-4.949E  01 

.60 

6*0005 

05 

7 

.PP1E-C7 

5.728E-03 

.5.127E  01 

.70 

7.CCCE 

C5 

7 

.  635E-07 

2.580E-03 

-5.283E  01 

.PC 

? • 0  jCZ 

cr 

V  — 

7 

.416E-C7 

i.940E-03 

-5.423E  01 

.?c 

'=*?PQE 

05 

7 

.1S9E-C7 

1.674t-03 

-5.551E  01 

l.OC 

1 .0005 

C6 

6 

.S78E-C7 

1.46CE-03 

-5.670E  01 

l.ic 

1*1 oce 

C6 

6 

.755E-C7 

1.285E-03 

-5.781E  01 

1.3C 

■l.^'OOE 

C6 

6 

.535E-C7 

i.l40E-03 

-5.885E  01 

1  >30 

1.3005 

•CA 

6 

.31£E-C7 

1.017E-03 

-5.984E  01 

1.40 

l.'^OCE 

06 

6 

.iC6E-C7 

9.159E-b4 

-6.078E  01 

l.BO 

i.rcce 

oe 

.PC1E-C7 

8.534E-04 

-6.168E  01 

1-60 

l>i:CE 

C6 

c; 

.7C3E-C7 

7.461E-04 

-6.253E  01 

1.7C 

1.7C;0E 

C6 

c; 

w 

.513E-C7  . 

6.788E-04 

-6.335E  01 

1  .so 

1  .="0E 

C6 

c 

.331E-C7 

6.198E-04 

-6.414E  01 

1-.90 

l.^CCE 

C6 

5 

.156E-07 

5.6S0E-04 

-6.490E  01 

?.cc 

•r.rccE 

06 

4 

.990E-C7 

5.222E-04 

or 

p»l- 

.f-^^CE 

C6 

4 

.SS3E-C7 

4.817E-04 

-6.633E  01 

S » 50 

? . ECCE 

06' 

4 

.6P3E-C7 

4.456E-04 

-6.7C1E  01 

5.3C- 

?.7CCE 

C6 

4 

.545E-C7 

4.133E-04 

-6.766E  01 

5.40 . 

?.4CCE 

C6 

. . 4 

.406E-C7 

3.844E-04 

-6.829E  01 

5.5C 

2.5CCE 

C6 

4 

.2P1E-C7 

3.584E-04 

-6.890E  01 

g.6C 

3.iCCE 

C6 

4 

.161E-07 

3.349E-04 

-6.949E  ^ 

2.7C 

P.7CCE 

06 

4 

.047E-07 

3.137E-04 

-7.006E  01 

.... 

■  'b.'scce' 

06 

•  .  3 

.939E-07 

2.944E-04 

-7.061E  01 

?.90 

2.S00E 

C6 

3 

.837E-07 

2.769E-04 

-7.114E  01 

3.00“ 

3.C:0GE 

C6 

3 

.74bE-C7 

?.'609E-04 

-7.166E  01 

3.10 

3.10CE 

C6 

3 

.645E-07 

2.461E-04 

-7.216E  01 

3.?0 

3.200E 

06 

3 

.5S6E-07 

2.326E-04 

.7.266E  01 

3.30 

3.3C0E 

C6 

3 

.470E-07 

2.201E-04 

-7.313E  01 

3.40 

3.4QCE 

06 

3 

.389E-b7 

2.b86E-b4 

■  -7«36bE  01 

3.50 

3.500E 

C6 

3 

.312E-07 

1.980E-04 

-7.405E  Oi 

3.60 

3.6C0E 

06 

3 

.S37E-07 

1.882E-04 

■;7.449E  01 

3.70 

3.700E 

06 

3 

.167E-07 

1.791E-04 

-7.492E  01 

3.80 

3.80CE 

C6 

3 

1.707E-04 

3.90 

3.90CE 

06 

3 

.C35E-07 

1.629E-04 

-7.575E  01 

4.00 

4.000E 

06 

5 

.973E-07 

1.556E-04 

-7.615E  01 

B.  Example  2 

Repeat  the  above  example,  but  quiet  the  rudder  discharge.  (This 
might  be  done  to  investigate  the  effects  of  adding  p-static  dischargers 
to  the  rudder  assembly  of  the  aircraft.) 


1.  Input  Deck 


The  input  deck  required  to  evaluate  this  problem  is  as  follows 


1 


3«TAILCAP 
,_A5i . 


+o.aie;-o3 

♦0.23E-53 

+  ^■•35^  -C 1 

C  •'•mz 

tailcap 

+0.35E-03 

+:.35E-01 

1 

tajlcap 

♦0.2CE-0.? 

+0»56E"?3 

♦0.351-01 

2  ^’HZ 

tailcap 

+0.3CE-02 

+0.16E-0? 

*3*35E-C1 

3  VH2 

tailcap 

+0 • S0E*0? 

+0.21E-0? 

+3*3?E-C1 

A  ''HZ 

■tailcap 

+0.27E-02 

+  C.11E-''? 

*c.37r-oi 

5  VMZ 

tailcap 

+0.27E-02 

+Ct75E-03 

♦o.^on-ci 

6  ''HZ 

tailcap 

+0.32E-02 

+0.1CE-0? 

+  :'.3?r-ci 

7  VHZ 

tailcap 

+6.A3E-0? 

+0.17E-C? 

+3.38E-C1  ^ 

F  “HZ 

tailcap 

♦0.70E-02 

+0>1CE*0? 

♦  C • 35E  *0 1 

?  VHZ 

tailcap 

+0.10E-01 

*C#36E-'0l 

10“HZ 

TAIlCA'p 

+C.13E-C1 

'►r*40E-Cl 

U“hZ 

tailcap 

♦6.13E-01 

+:.E1T-C1 

12“hZ 

tailcap 

+0.12E-01 

•*»0«9CE-73 

+:.^7E-C1 

13“hZ 

tailcap 

'  '+0.102-01 

+  0#1CE*C.2 

14“.HZ 

tailcap 

It 

2 

KC-135 

1.00  600*0  ?C.C 

0 

0.10  4.0C  O.IC 

8.6  t  1*0 _ 

1»CIPRU3  CL‘’'J0 
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2.  Output  Deck 
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to  (Ai  fi)  caM<4>  o)  ro  i\>  j\)  ^oIiO  v  'vj  >\»  i\) 


S3 1 

P-STATIC  MeSEL  CC8NTD] 

■"  CALCU 

LAT^rO  C^U.RGT^.3 

CURRENT  IS  i.oocr-03 

AMPS 

9ILITY  IS  .CC2 

0  TWAT  The  CWAR3ING 

CURRENT 

WILL  = 

E  OPrATO-  T^AA: 

l.CCCr-03  Ani=>S 

OHeRT-CI»CvlT 

equivalent 

NeiSE  FIELD 

1^7] 

CA^=^S1 

CVSLTS/M3 

.10 

l.OOOE  Or 

1 ,B49E-C8 

3.S70E.^04 

•  ?0 

?;.;ccE  C5 

i.p34F»C8' 

1.920E-04 

7.0COS  C5 

1.612E-C8 

1.264E-04 

•  40 

A.OOCE  C- 

1.784E-0S 

9.334E-Q5 

O.OOCE  Or 

1.75CE-08 

7.325E-0f 

•  f'O 

O-.COOE  05 

1.711E-C8 

5.97CE-05 

•  70 

7.: OOF  05 

1.67CE-CR 

4.992E-05 

.80'  '■ 

. --c.y5e---er- 

1.6c5E-C8 

4.252E-65 

'••O  jOt  05 

l.EeOr-C£ 

3.674E-05 

i.cc 

i.5:^cce  06 

1 .534E:*CS 

3.210E-05 

1.10 

1.1 :cE 

1.510t-C8 

2.873E-0? 

l.'OCE  C6 

1.489^-08 

2.597E-05 

1^3C 

l.'OCE  C6 

1 .ATEr-CS 

2.36SE-05 

i  *^0 

. “C^T" 

1  .‘‘975-08 

2.178E-05 

1 

1*"‘CCE  C6 

1 »  445“-08 

2.016E-05 

1*60 

1-^=:0E  06 

1.4395-08 

1.873E-0E 

1^70 

i«7:GE  C6 

1 .4EfeE-c« 

1.758E-05 

1 

.  l.^OCE  06 

1.4P2E-C8 

1.693E-05 

1»9C 

06 

1.417E-CS 

1.561E-05 

.^•00 

".:dce  g6 

1.480E-6E' 

?*10 

E.l'CE  06 

1.79CF.-C8 

1.784E-05 

'.r:cE  C6 

2.215E-C8 

2.108F--05 

?*30 

5.0- -^CE  C6 

2.65CE-C8 

2.411E-05 

?*40 

^.‘*GOE  06 

3 . 073E-C8 

2.6S0E-05 

^•50 

5. 5 COE  C6 

3.4P1E-C8 

2.914E-C5 

5.6::CE  C6 

"  3.?7'6E-C8 

. S.lTSE-dE 

?.7G 

S.TGCE  C6 

4.240E-08 

3.286E-0E 

^••«0 

5.^0CE  C6 

4,59CE-08 

3..431E-05 

2*90 

E.r-oci  C6 

4.9235-08 

3.553E-05 

3.00 

3.C0CE  06 

3.238E-08 

3.654E-05 

3tlO 

?.10CE  C6 

5.3922-08 

3.640E-0E 

^.“COE  06 

5.538E-C8 

3'7g22E-dr 

3*30 

3.300E  06 

5.677E-08 

3.601E-0E 

3*40 

■3. 4 COE '06 '■ 

. "  5.  MCE -08 

3.576E-05 

3.50 

■ 

3.5CCE  C6 

5.936E-08 

3.550E-05 

3.6CCE  Q6 

6.056E-d8 

3.521E-05 

3.70 

3.7.CCE  06 

6.170E-C8 

^3.490E-05 

ECU  I VALENT 

NeiSE  FIELD 

tDSV/M] 

-6.8a3E  01 
-7.A32t  01“ 
-7.795E  01 
-8.05SE  01 
-5.269E  01 
-8t447E  01 
-8.6C2E  01 
-8.74 IE' 01 
-8.863E  01 
-S.985E  01 
-9.082E  01 
-9.169E  01 
-9.249E  01 
-9.322E  01 
-9.389E  01 
-9.451E  01 
-9.508E  01 
-9.562E  01 
-9.611E  01 
-^.6g8E  01 
-9.495E  01 
-9.350E  01 
-9.234E  01 
-9.142E  01 
-9.069E  01 
■-9.diiE'oi 
-8.965E  01 
-8.927E  01 
-8.897E  01 
-8.873E  01 
-S.876E  01 

. “-fi'gi'orur" 

-8.886E  01 
-8.891£'01 
-8.898E  01 
■  -SisoSE  01  ' 
-8.913E  01 


C.  Example  3 

Calculate  the  equivalent  noise  field  induced  in  a  belly -mounted 

antenna  on  an  F-4  aircraft.  Assume  that  the  antenna  has  an  induction 
2 

area  of  8.6  m  ,  and  that  the  receiver  has  a  bandwidth  of  1  kHz.  Further 
assume  that  the  aircraft  is  flying  at  a  speed  of  600  mph  at  20  kft  through 
stratocumulus  cloud.  Allow  all  extremities  of  the  aircraft  to  discharge 
and  evaluate  the  ENF  at  uniformly  spaced  frequencies  of  0.1  to  4.0  MHz 
with  a  Af  of  0.1  MHz.  (The  F-4  is  approximately  1/3  the  size  of  a  KC-135.) 

1.  Input  Deck 

The  input  deck  required  to  evaluate  this  problem  is  as  follows: 


1 - 

O-BELWY 

15/ 


+0.14E-03 

+0.30E-03 

+0.90E-C5 

♦0.15E-03 

+:.11E-C3 

+C.20E-C3 

♦0*r7P*03 

+C).16E-0? 

+  0 1  ??5e*05 

+  •  s  0  •  b  3 

+C.1CE-0? 

+0.17E-02 

+c,4or-C3 

+0.3CE-03 

+C.12E-C3  . 

+0.50E-C3 

+0t55E»03 

*  - .23E-03 

+C.85E-03 

+O.11E-02 

+0.401-03 

40.17E-02 

+0.27E-C? 

+C'»  icr  -  c? 

•  ?  9  E  -  'S  ^ 

•  1  W  7-  •  4.  ^ 

+0.22E-02 

+0 • 39^-03 

16E-C3 

+C.15E-02 

+  :’ •  lOE-C? 

+C»1SE-C2 

+:.70E-:3 

■'+0.i9E-C? 

+b*B0E"O? 

+  C‘«63E"03 

+0.20E-02 

•^0«46E-C? 

+:.6cr-c? 

T7 

1 


F-4 'fighter 

0*33  6C0.0  SC.C 

0 . 

0«10  4»00  0.10 

S.6  /  1.0 

2«STt?AT5  CU 


0  3EULY 

I  ^HZ  3E1.LY 

C  mmz  belly 
'■y^Hz  ??lly 

4  !-'.HZ.  belly 

5  belly 

6  MHZ  belly 

7  MHZ  belly 

bmhz  belly 
■y-WZ’-BECXT' 
10  MhZ  belly 

II  mhZ  belly 
13  MHZ  belly 

13  mhZ  belly 

14  MHZ  BELLY 
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2 


Output  Deck 


The  program  output  is  as  follows: 


S^I  STATIC  ^Lg:CTRIClTY  MSDEl 


P-STATIC  ^1'^DEL  FVA^JATED  P£R  A  F-4  FI3HTEP  AIRC^^AFT 
WITH  The  PECEIVINj  A-JTENNA  located  at  the  3ELLY 


_ scale. SIZE 

S-E^D 

altitude 

CKFT] 

CLaUD  TYPE- 

. #33  ■ 

600  *0 

20*0 

strats  CU 

"'““START  RREQ. 

STop  FREG* 

DEuTA-cr 

C'^HZI 

C'li^ZD 

rYHzi 

.10 

4.CC 

.10 

RECEIVER 

^eisE 

pandwidth 

antenna 

INDUCT  I aN 
a^ea 

ExhIT 

■  j 
i 

1 

t 

i 

!  • 

O 

O 

^  0 

A'LL  DISCHARGfTS 

PEP^^ITTEO 
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_ S£l_ 

P-STATir  MftOn _ CCaNTDl  _ 

.  ■  -..THE  ...CALCJ 

LATPC  C'-<4RGIN'5  CUREEM  IS  6»6.00E- 

04  .AMP5, 

jwr  OR03A 

3ILITY  IS  .0061 

THAT  THE  C^ARGIN 

3  CURRENT 

WILL  B 

E  Creates;  than  6 

.6C0E-04  AHPS 

FCrojEKCY 

SneRT-CIPCylT 

equivalent 

equivalent 

CU^^'EMT 

NOISE  FIELD 

N9ISE  FIELD 

C'^WZI 

c-z] 

CANPSl 

tveLTS/M] 

CDBV/HJ 

. .‘IQ^ 

■  lOOCE  05 

1.171E-07 

2.45bE-03 

-5.221E  01 

•  ?d 

"•COCE  05 

1.166E-C7 

1.221E-03 

-5.826E  Oi 

•  30 

9*Z-'^0Z  05 

1.157E-07 

8.073E-04 

-6.1S5E  01 

•  40 

4*r..0CE  Or 

l.lA3r-C7 

5.983E-04 

-6.445E  01 

•  SC 

S*"noE  Or 

1.126E-07 

4.712E-04 

-6.652E  01 

•  60 

■f.rC'CE  05 

1.1C4E-07 

,  3.853E-04 

-6.827E  01 

.70 

T-'OCE  0*=^ 

1.C51P-C7 

3.E3EE-04 

-6.980E  01 

.•?0 

^•COOE  05 

1.C55E-C7 

2.766e-64 

-7.117E  01 

•  '^0 

■’♦COCE  C5 

1.D3SE-07 

C.391E-04 

-7.242E  01 

1  •oc 

l-COCE  06 

1.0GCE-C7 

J?.C94E-04 

-7.357E  01 

1»10 

l*lr.OE  06 

9.722E-C8 

1.850E-04 

-7.464E  01 

1  *20 

l.IOC-E  C6 

9.441F-08 

1.647E-04 

-7.565E  01 

.  It30 

i»~“;e  c6 

9.1642-08 

1.475E-04 

-7,661E  01 

i.itc 

i.^>nE  C6 

B.893G-GS 

i.336E-04 

-7,751E  01 

l.SG 

1  •  5  r  '  i  ^ 

?. 6395-08 

1.204E-04 

-7.837E  01 

1  *60 

l.^'CE  06 

8.374:-0S 

1.095E-04 

-7.919E  01 

1  .70 

1*700E  C6 

«.12S“-08 

1.001E-C4 

-7.998E  01 

!  .BC 

1*':CE  C6 

7.893E-CS 

9.178E-C5 

-S.073E  01 

1.90 

l.?OCE  C6 

7.667E-08 

8.446E-05 

-8.145E  01 

5*00' 

E.'OCE  C6 

7.AE2F-og 

7.798E-CI5 

-8.215E  01 

?.10 

E.IOOE  C6 

7.,?A6E-08 

7.222E-05 

-8.281E  01 

?.2C 

3*5C'CE  06 

7.C50E-08 

6.707E-OF 

-3.345E  01 

?.30 

?*2CG£  C6 

6.863E-0S 

6.246E-05 

-8.407E  01 

?.40 

2.!*GCE  06 

6.685E-08 

3.S30E-0S 

-8.467E  01 

?.50 

F.50CE  06 

6.E16E-08 

5.456E-05 

-8.525E  01 

2*60 

'^OCE  C6 

6.556E-08 

5.116E-b5 

4|,-§8ie  W" 

?.7C 

?*7CCE  C6 

6.2C2E-08 

4.808E-0E 

-S.634E  01 

?tSO 

3*?0CE  C6 

6.0E7E-08 

4.527E-05 

-S.687E  01 

2*90 

3.9CCE  C6 

5.918E-C8 

4.271E-05 

-8.737E  01 

2*00 

3.00CE  C6 

5.786E-C8 

4.037E-0F 

-8.786E  01 

3»10 

3.1C0E  C6 

5.678E-08 

3.834E-05 

-9.831E  01 

3*20 

5.?0CE  C6 

S.583E-08 

3“.i52E-05 

3. 30 

3.3C0E  C6 

5.4‘92F-08 

3.483E-05 

-8.914E  01 

3ir4'o . 

3.4CCE  C6 

5.4b5E-08 

3.327E-05 

-8.954E  01 

3 150 

3.50CE  06 

_ 5_*3?..1.^.-Q8.._ 

3.182E-05 

-8.993E  01 

3 '.60 

3.66cE  06 

5 .2412-08 

3.b47E-b'5 

•9.031E  01 

3.70 

3.70CE  06 

5.164E-08 

2.921E-05 

-9.067E  01 

3.SCOE  06 

5.05lfe-08 

2.S04E-05 

3.90 

3.30CE  06 

5.030E-08 

2.694E-05 

-9.137E  01 

4.00 

4.C00E  06 

4.952E-08 

2. 59 IE -05 . 

...  Qj- 
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D.  Example  4 


Repeat  Example's,  except  assume  that  the  aircraft  is  now  flying 
through  cirrus  cloud. 


1.  Input  Deck 


The  input  deck  required  to  evaluate  this  problem  is  as  follows: 


—1 - ; - ^ - 

._cL-mLy. 

15* 

+0.14E-03  +0.?CE-03 

+-'?«90E.*05 

0  ^HZ  3ELLY 

♦0.15E-03  +0.22E-C3 

+0.111-03 

1  hhz  belly 

♦C.20E-03  +O.E7E-03 

+o.is-:-03 

2  “^HZ  belly 

♦0.i6E-Og  *0.522-03 

3  VHZ  belly 

+0.1CE-02  4.0.17E-02 

+:.*0E-C3 

4  HMZ  belly 

+Cil0E-03  +0.802-03 

+0.12E-C3 

5  -IMZ  BELLY 

♦0.5CE-03  +0.5EE-03 

+  j  #  •’03 

6  YHZ  belly 

+C.85E-03  +0.11E-02 

+  0.40'"-03 

7  MMZ  B2LLY 

♦0.17E-02  +0.27E-02 

+0. ICE -02 

8MHZ  Br-LLY 

♦0.2^E-02  +Q.29E-02 

+  0.l5t--C2 

9  ’^HZ  BELLY 

+0.22E-02  +0.29E-02 

+0. 16E -C2 

10  HHZ  BELLY 

♦C.15E-02  +d.42E-02 

+0.10E-02 

11  MHZ  BELLY 

+0.18E-02  +0.6EE-02 

+0.70E-C3 

12  MHZ  §ELLY 

+0.19E-02  +0.5CE-0? 

+0.62E-C3 

13  MHZ  BELLY 

♦0.20E-02  +G.46E-02 

+  0 . 60c.  -03 

14  MHZ  BELLY 

li 

1 

F-4  FIGHTER 

0.33  600.0  20.0 

0 

C.IO  4.00  0.10 

8.6  t  1.0 

1»CIPRUS  .CI,6U?. 
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2.  Output  Deck 

The  program  output  is  as  follows: 

S'?!  STATIC  electricity  MeCEL 


P-STATIC  .evaluated  FeR  A  F-4  FIGHTER  AIRCRAFT 

WITM  thc  RECEIVIVG  A^'TEN^A  LeCATFC  AT  THE  BELLY 

"altitude  "CLSuiTTYPE 

CKFT] 

20tC  CIRRUS  CLeUD 


SCALE  SIZE 


SPEED 

[YPH] 


.53 


6  jO  •  C 


STA5T  FREQ. 
[MHZ  3 


ST5P  FREQ. 
[MHZ] 


.10 


4. CO 


CELTA-F 

[MHZ] 

.10 


RECEIVER 

N9ISE 

BANDWIDTH 

CKHZ3 


antenna 
INDUCT  I SN 
area 

CH*»23 


.L»_00  S.60 


ALL  DISCHARGES  PERMITTED 
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SRT  o-STtTir:  ’^ssEi.  ccasiTn] 


TMr  CALCULATE?  CHA(53 

l\3 

CyRRE’^T  IS  3.3C0E-0A 

AMPS 

THE  pRSBABILH'^  is 

.CC61  THAT  THE  CHARGING 

CURRENT 

WILL  3E 

GprATOP  T 

-:A\ 

3.3C0E-CA  AH?S 

TQr  r:uENcv 

\CY 

SH?RT-CI5CyIT 

equivalent 

equivalent 

CyRRE^:T 

NSiSE  FIELD 

N8ISE  field 

r-izi 

[AHBS] 

CveLTS/M] 

[DBV/M3 

•  ID 

l.OCCt 

rC. 

8.a77E-08 

i.732E-03 

-5.522E 

01 

V20 

r  ..:g:e 

S.EA7E-0S 

3*63lE-04 

-6.127E 

oi 

#3C 

O'z 

S.IBSE-CS 

5.7C9E-04 

-6.486E 

01 

A  •  I OCE 

r  cr 

V  -- 

s.C?5E-C.?.  ' 

A.230E-04 

“6.7A6E 

01 

•sc 

S  •  C  GCE 

OS 

7,059?. 03 

3.332E-04 

-6.953E 

01 

.60 

4.:.;nE 

7.P.1CE-C3 

2*724F-04 

-7.12SE 

Cl 

•  70 

cs 

7.6A2E-C3 

2.285E-04 

-7.281E 

01 

^  .?o 

SOOGE 

r  5; 

V  — 

7,46iE-C8 

i ♦ 952Et04 

-7,4i«E 

01 

•  9C 

05 

7,27CE-C8 

l*691E-04 

-7.543E 

01 

!*CD 

1 • OOCE 

'"6 

7.C73E-C8 

1*480E-C4 

-7.658E 

01 

1.10 

1  .  lOOE 

06 

6,874E-C8 

1.303E-04 

-7.765E 

01 

1 .30 

1  •  ?  Ov'^E 

6  *676 E “08 

1.164E-04 

-7.866E 

01 

1.30 

_ l.!liS^ 

•'^6 

. 0.8 

1*043E“04 

-7.962E 

01 

1.4C 

i  .^OCE 

C6 

6 • £88E”C8 

3.401E-0E 

-8.C52E 

01 

l.SO 

1*SjOE 

C6 

6.1C2E-CS 

S*S14E-0E 

-8.138E 

01 

1.60 

1 -rOCE 

06 

5.921E-C8 

7.745E-0H 

-8.220E 

01 

1.70 

1 .7QQr 

06 

S.748E-C8 

7.C77E-05 

-8.299E 

01 

1.^0 

1*'00E 

06 

5.5S1F-C8 

6.450E-05 

-8.374E 

01 

l.OG 

1  .t^^OE 

06 

?.4?2E“C8 

5*972E-CF 

«S*446E 

01 

?'.oo 

?*:doe 

“5i" 

5'.269E-C8 

■■  5t514E-0  = 

-8.516E 

01 

■  3*10 

?*ncE 

06 

S.124E-C8  , 

5.107E-0E 

-8.5S2E 

Cl 

?.?o 

2»'CCE 

06 

4.SS5E-C8 

4*7432.05 

“8.646E 

01 

?.30 

2.3CCE 

06 

4 .8S3E-08 

4*416E-0S 

•8.708E 

01 

?.40 

?  •  AOOE 

06 

4.727E-08  - 

4.123E-0S 

-8.76SE 

01 

?.50 

2  .SOOE 

C6 

4 .6C8E-C8 

3.858E-0E 

-8.826E 

01 

?>.60 

2.6:cE 

06  ' 

4.494E-08 

3.6iSE-6E . 

■  ■ -8.882E' 

oi 

3. 70 

06 

4.356E-C8 

3.400E-05 

•8.935E 

01 

?.80 

i*SCCE 

06 

4.283E-C8 

3.201E-0S 

-8.988E 

01 

2.90 

2.P0GE 

C6 

4.185E-0S 

3.C2CE-05 

-9,038E 

01 

3.00 

3.DOCE 

06 

A.091E-0S 

2*85AE-0E 

•■9.087E 

01 

3.10 

3.1G0E 

06 

A.015E-0S 

2.711E-0E 

-9,132E 

01 

3.20 

3.200E 

C6 

3.548E-08 

2.E82E-05 

-9.174E 

01 

3.30 

._.  ?.?3CCE,. 

06 

3.SS3E-08 

2.A63E-05 

-9.215E 

01 

3.40 . 

3»**00E 

C6 

3.822E-08 

2.353E-05 

-9.255E 

01 

3.50 

3.5C0E 

C6 

3.763E-08 

2.250E-95 

-9.294E 

01 

3.60 

3.60CE 

C6 

3.7C6E-08 

2.155E-05 

-9.332E 

01 

3.70 

3.70CE 

C6 

3.652E-C8 

2.066E-05 

-9.368E 

01 

?T?0 

3.5C0E 

C6 

3.600E-08 

1.983E-05 

-$.A04E 

"or 

3.90 

3  •  900E 

06 

3.550E-08 

1.905E-05 

-9.438E 

01 

4.06 

4.C00E 

06 

3.562E-C8'^ 

i .832E-05 

-9.472E 

01 
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E.  Example  5 


Using  the  streamering  model,  evaluate  the  ENF  induced  in  an  antenna 

mounted  near  the  radome  of  a  B-47  bomber  due  to  cirrus-cloud-caused 

p-static  charging.  Assume  that  the  antenna  is  0.04  m  aft  of  the  front 

2 

of  the  radome,  and  that  the  antenna  has  an  induction  area  of  0.01  m  . 
Assume  that  the  minimum  characteristic  dimension  of  the  radome  is 
0.24  m  and  that  the  ratio  of  the  dielectric  frontal  area  to  the  total 
aircraft  frontal  area  is  0.01.  Further  assume  that  the  size  of  the 
B-47  is  0.89  times  the  size  of  a  KC-135,  and  that  the  B-47  is  flying 
at  600  mph  at  20,000  feet  through  cirrus  cloud. 

Evaluate  the  ENF  at  nonuniformly  spaced  frequencies  of  1.13,  2.16, 
4.35,  8.62,  and  10.7  MHz  for  a  receiver  noise  bandwidth  of  1.0  kHz. 

Deck 

The  input  deck  required  to  evaluate  this  problem  is  as  follows: 


”2 - 

Nk  RA0<9vE 
b-47  eSMPER 
C.S9  600.0  ?C.O 

1 

5 _ 

-►1.13E  +  CC 
^8#16E^C0 
♦  4.35E«**00 
♦St62E^00 
• 07E+01 
Otic  1.00 
1*CIRRUS  • 
2«RADeME 

C«C4  0t?4  0#30 

C.Ol 
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2.  Output  Deck 


The  program  output  is  as  follows : 


S=^I  static  electricity  M8DEL 


P-STATIC  l^QDEL  EVALUATED  PER  B-47  BSMBCR  AIRCRAFT 
WITH  THE  RECEIvTNj  ANTEN!\!A  LSCATEd  AT  THP  NR  HAOSH 


FBR  STREAMERING  CCCuRING  EN  THE  RAD?^^E 
AND  the  antenna  , OA  HETERS  '  APT  BF  THE  PR8NT  BP  THE  RAD5ME 

AND  A  HINIHUM  characteristic  DIMENSI5N  0P  .24  METERS  8P  THE  DIELECTRIC  RAOOME 
AND  A  FUSELAGE  DIAMETER 'SP  '  .30  METERS 

AND  A  DIELECTRIC  AREA  TS  A/C  FRSNTAL  AREA  RATIE  8P  »01 _ . _ 


SCALE  SIZE 

S^EED 

[MPH] 

altitude 

C<PT3 

CLOUD  TYPE 

.SS 

60C*0 

2C.0 

CIRRUS 

start  preg. 

—  "jrqrQ"  ' 

delta-p 

[MHZ] 

C^^HZ3 

tMHZ3 

1.13  _  _  .4.0 ‘70  Nes-yMPeRM 


RECEIVER _ ■  antenna 

NOISE  INDUCTION 

BANDWIDTH _  area 

CKHZr  ■  ■  ■  tM**23 

iToo'  . ■■■'.10 
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THr  CALCULATrP  CHARGING  CURREKT  IS  8 •900^-04  AMPS  . 

THE  PPOBABILITY  IS  «C02g  IHAT  THE  CHARGING  CURRENT 
WILL  BE  GREATER  Ti-iAN  A.goOE-OA  AMPS 


THE  CALCULATED  STREAMER  I NG  CURRENT  IS  8*90E-C6  AMPS 


frequency 

FRE3UENCY 

shrrt-circuit 

CURRENT 

EQUIVALENT 
NSISE  FIELD 

EQUIVALENT 

NOISE  FIELD 

CMHZ3 

CHZJ 

CAMPS] 

CV0LTS/M] 

CDBV/M] 

TiTs 

T.isCE  'c6 

6.29SE-10 

1.CC3E-04 

-7.996E  01 

_ g«16 

_ _ £_*J^.G£..C6 . . 

,  ..  2.292E-1C 

1.910E-05 

•9.436E  01 

4.35 

4..35CE  C6 

6.883E-11 

C.8.4SE-06 

-1.109E  02 

8t62 

8.62GE  06 

1.913E-11 

..  3.995E-07 

10.70 

1.C7CE  07 

1.P60E-11 

2.119E-07 

-1.335E  02 
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Appendix 


PST  AT  PROGRAM  LISTING 
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*.°SJAX* _ eCT  1974  VE^3IfiN  D»*g  SRI>  MENL9  PAg,l<:j_CA[_« 


- PSTA.T_C9MPJT!;5  equivalent  \-9ISE  g'lELDS  GENERATFn  IN  AN  P' 

A1RC9APT  ANTENNA  DjE:  T3  ELECTROSTATIC  DISCHARSES  8CCURING  FR9M  ThE - P 

airfoil  EXTprvtTIES.  oc; . .  ,  p, 

PSTAT  CS'^PiJTrs  TrJE  E3'JlVAi.ENT  N8iSE. FIELDS  INDUCED  IN  AN  AIRCRAFT  P; 
ANTENNA  DyE  TS  STREayr:;!^^  DISCHARGES  ?N  DIELECTRIC  CANBPY  OR  P' 

radohe  surfaces*  pi 

_ THE.  USER  can  select  EIThER  ngpE  sp  PRSGRAH  EXECUTION  BY  AN  Pi 

apfrgbtate  data  card*  pj 

PRESENT  (1974)  CojolINj  DATA  (DATA  DESCRIBING  THE  ELECTROMAGNETIC  Pi 

Coupling  between  A^  airfoil  ti=  and  an  antennai  snly  exists  for  P' 

belly-  and  TAILCAB-,<D!jNiTED  ANTENNAS  AND  DISCHARGE  LOCATIONS  AT  THE  Pi 
and  elevator  tips*  (STHER  possible  DISCHARGE  LOCATIONS  PI 
ARE  UMH=«RT4nt  pg;;  ppisONS  DESCRIBED  IN  THE  FINAL  R'EPORtT*  Pi 

THE  PROGRAM  is  GENERAL  I ZED^  SO  THAT  AS  ADDITIONAL  COUPLING  DATA  pi 
becomes  AVMLA=LE>  .IT  may  be  INCORPORATED  INTO  THE  PROGRAM*  THr  P' 
additional  DATA  *'AY  DE  AN  EXTENSION  OF  THE  FREQUENCY  RANGE  OF  THE  Pi 
-g.Xl.STTN3  DATA  (IV  l-»p.7  intervals*  U?  T9  100-MHZ)j  SR  COUPLING  PI 

DATA  (ASAr-',  IN  intervals.  UP  TO  ICO-mhZ)  FOV'antENNAS’ .  . 'P! 

located  IV  other  positions*  the  CSU®L1NG  DATA  USED  IN  PSTAT  IS  P! 
EXPEPIMENTAL  DAT.a  “^yTAlNED  pROh  :<C-135  SCALE  MODEL  AND  FLIGHT  TESTS.  PI 
AND  IS  READ  IKT^  THf  PROGRAM  pcjgy  CARDS*  p( 


extremity-to-tailcab  COU 
belly  (PUSELAGErCOyPLiN 
AOPROOIATE  TO  H!£  NEEDS* 


SINCE  THE  SPECTRy 
A  IDO-A'lHi  prEQUENCV 
EST.  AND  PSTAT  ORES 
'  Below  To'O-^wzv  “sh? 
PROGRAM  MgDIEiCATlS 
USERS  GUIDE  ROR  DlR' 
DUE  TO  THE  NATJR- 


M  OP  CORONA 

“range  is  aD 
ently ^limits 
u'lc  a  higher 
N  may  be  mad 
ectisns*"" 

■“  STREAMER 


PROGRAM  ARE  DEVOTED 
NOISE*  the  desired 
DATA  CARO  READ  INTO 
SECTION  ONE.  THE  CO 
2.  THE  STREAMER ING 


TO  the  calc 
SECTION  IS 
THE  PROGRAM 
RONA  SECTION 

section* 


ULD  SELECT  THE  DECK 


DISCHARGE  NOISE  FALLS  OFF  AS  1/F. 
equate  to  handle"  M.6§r'CASE§  OF  INTER- 
_THE  CALCULATION  TO  FREQUENCIES  AT  OR 
FREQUENCY  '^ANGE 'BE ‘^OES IRED.  A'  SIMPLE' 
E  TO  DO  SO.  AFTER  CSNSULT_ING_THE _ 

ING.  AND  THE  INPUT  REQUIREMENTS  FOR 

ELDS,  separate  SECTIONS  Of  THIS - 

ULATION  OF  streamer  NOISE  OR  CORONA 
SELECTED  BY  THE  USER  AS  THE  FIRST 
*  A  1  (ONE)  ON  INPUT  IMPLIES 

A  2  (TWO)  ON  INPUT  IMPLIES  SECTION 


PSTATOOl 

PSTAT002 

PSTAT003 

PSTAt004 

PSTAT005 

PSTATOOO 

PSTAT007 

PSTAT008 

PSTAT009 

PSTATOi() 

PSTATOll 

PSTAT012 

PSTAT013 

PSTATOlO 

PSTAT015 

PSTAT016 

PSTAT017 

PSTATOia 

PSTAT^B 

PSTATOa^- 

PSTAT021 

PSTAT022 

PSTAT023 

PSTAT024 

PSJAT025 

PSfATOaS 

PSTAT027 

PSTAT028 

PSTAT029 

PSTAT030 

PSTAT031 

pstTtW 

PSTAT033 

PSTAT034 

PSTAT035 

PSTAT636 

PSTAT037 

PsTAmo 

PSTAT039 

PSTATOOO 

PSTAT041 

PSfATd42' 

PSTAT043 


PSTAT045 

PSTAT046 

PStAT047 

PSTAT048 

PSTAT049 


PSTAT051 

PSTAT052 

PSTAT053 

PSTAT0S4 
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_c_ 

c 

e 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

"r 

c 

c 

c 

c 

c 

c" 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r 

c 

c 

c 

c 

c 

IT 

c 

c 

c 

c 

c 

T“ 


- - lE-Lf.*;!/  IS  3N  .(gP  NEAR)  ELEVAT9R  TTP 

IP  LA»?/  A\Tr\\A  IS  ?Ni  (3R  N'SA^)  WING  TIP 

.  l-A"3/  AVTtNNA  IS  !3Ni  (?R  MEAR)  j^ING  TIP 

LANT»14  character  AL»HANUvie:RIC  CESCRIPtlSN  OF  ANTENNA  L3CATieN 

IERR»  ERRSR  plAG--  3ET  =  1  IF  DATA  INPUT  E=?R8R  9CCURS 
EP3IL=  E'=SIL'*^--  =E'?''1TTIVITY  5F  FREE  SPACE  (FARADS/METER) 

C9U=LiNS  Coefficients  Te  be  read  (Ncaup  alss 

DEFINES  Th£  MAxif^U"'  FREGJE.NCY  +  1NHZ)~  ’  " - 

EST9jWST'>^RST‘?*  storage  ARRAY?  PSR  NCSyP  C8UPLIN3  C9EFFICIENTS 

=  RD''  ELEVATeRSj  *INGS/  RUDDER  T5.  SELECTED  ANTENNA 
v'^CATIDN 

NC9»  NC9JP  +1 

CYCLES  T£  3E  MADE  USING  THE  SAME  C8UPLINS 

DATA>  ?U'^  (i^fSSISLY)  VAFI?US  FARANET ERS" . . . .  . 

I8PP=C8R9NA  DISCHARGE  DUENCH  CeDE  (AIRF0IL(S)  P-STATIC  PR8TECTED) 

*  I --ALL  discharges  PER^'ITTED 
».?--RUD^ER  ;I3C“ARGE  GUIETED  ?Y  4C  Di? 

=  3--NING  TI-‘’S  DISC=hAR3E  GUIFTED  3Y  40  DB 
. . 3UIETE0  8Y_40  OB 

*  E-“»JDDi--'5  A'D  ••••ING  "TIP?  DISThARG^'  G'^IETED' 3Y' ifO  'b 

*  6--RJDDE.R  A\D  ELEVATSP  t  ibc.  d  i  SchaRGES  DUIETEO  BY  40  D3 
»  ’--FLEVATeR  a\-D  XING  TIPS  DISC'^ARGES  S'UIETEO  BY'  40  DB 

IT«  6  WORD  al°'-A\.j«ERIC  CESCRIPTI8N  5F  AIRCRAFT 
XN«  aircraft  scale  size  (RELATIVE  T8  A  i<C-l3F) 

SPps  aircrapt  speed  (  in  ’^'iLES/hBjR) 

A'LT-  aIrCRAPT  ALTItlD'E  ■  ("In  KILSEe’ET')"' . . . 

Mei^cF.  PRrQ.jrvDY  SELECT  ^^CE  (.EG.  /  hfaNS  uMFSRM  FREQUENCY 

INTPPVA'.Si  .\P.  D  MEANS  USER  SELECTED  FREQUENCIES#  UP  T8 

FSTRTb  S'i'ART  prEGJENCY  (In  )  IF  yfDEF  .EQ«  C 
FSTP-STe^’  PPPOUENCY  (IN  IF  y^PEF  .ED»  0 

“DEL*  DELTA  r=EGJE\CY  (IN  MHZ3  IF  *'8DEF  .pqT  C . '  *  '  . . 

NPRr  number  <»p  FF-sUENCIES  T?  BE  EVALUATED  IF  MBDEF  ,ne,  c 

fppgj«  array  t3  Contain  user  selected  frequencies  if  “aoEF  .ne.  o 

Aant*  antenna  InDjCTI9N  Area  (IN  SQUARE  METERS) 

BNDW*  RECEIVER  NBicr  BANDvVlDtk  (IN  kMZ  ) 

(l«:iRRUS#  2.STRATS  CUMULUS#  4.FRSNTAL  SNSW) 

IC»  7  W5RD  alphanumeric  DESCRIPTIPN  "p  CL^UO  TYPE  (SEETCLST - 

CLSU*  FlBATInG-ocINT  ICL® 

SPDFA»  SPEED  FACtrR—  CHARGING  CURPE'^T  IS  RELATED  T8  AIRCRAFT 
„  . ,  S=EED  through  THIS  FUNCTI9N 

CH3C*  calculated  CHARGING  CURRENT  ( ^DISCHARGING  CURRENT)  (IN  AMPS) 
PR9B«  CALCULATED  PROBABILITY  «>F  CHARGING  »GT»  CHOC 
E#w#R«  Working  sisra^e  arrays  FbR  E:L?VAty!?#wlN6#  AN&  ftyDbER 


PSTATQ55 

PSTAT056 

PSTAT057 

PSTAT058 

PSTAT059 

PSTAT060 

PSTAT061 

PSfATO?F' 

PSTAT063 

PSTAT064 

PSTAT065 

PSTAT066 

PSTAT067 

"psTatCSs  ■ 

PSTAT069 

PSTAT070 

PSTAT071 

PSTAT072 

PSTAT073 

PSfAT074 

PSTAT075 

PSTAT076 

PSTAT077 

PSTAT078 

PSTAT079 

PSTAtOSO 

PSTAT081 

PSTAT082 

PSTAT083 

PSTAT084 

PSTAT085 

PSTAt1586'  ■ 

PSTAT087 

PSTAT088 

PSTAT089 

PSTAT090 

PSTAT091 

TgTiRT5?r- 

PSTAT093 

PSTAT(j94 

PSTAT095 

PStAT096 

PSTAT097 


C 

"C 

C 


C 

C 

T 


,.?9:^?yj1_H-EFRlCIENTS_lM8DIFIED  T8  .ACC8UNT  F8R  ANTENNA 

rNDUCTlON  ARE'A)  .  . . . . .  . . 

DISTRIBUTION  ^F- DISCHARGE  CURRENT  OVER  VARIOUS 
aircraft  EXTREMITIES 

D2R#02E#D2W«  DISCHARGE  CURRENT  SPECTRUM  NORMALIZERS 

xcou"  maximum  fresuenCV  5f  Coupling  data - - - 


■PBTAT098 

PSTAT_0?9 

PSTAtlOO” 

^TATlOl 


PSTAT102 

PSTAT103 


C 

C 


F»  FREQUENCY  CURRENTLY  BEING  EVALUATED 

LF»  COUNTER  FOR  FREQJ  - 

EX<  PRESSURE  COEFFICIENT  (P ( TORR ) »760»EX ) 


TsTATlOA 

PSTAT105 


ALPHA*  CORONA  PULSE  DECAY  TJME  CONSTANT 


PSTAT106 

PSTAT107 


PSTAT108 
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XNJ»  CSR^NA  “ULSE  RE-ETITIQN  raTE  - - 

. . T£ST_»_frequency  scaled  re  aircraft  _scale  size 

e^EGAs  RADIAV  ,p'RE3UENCY  . 

.......  PREL»^ELA.TIvr  p'jlSE  SPECTRu^*  AyPj. I TUDE 

b?MRiD3MEiiDe‘';/is  i  JSauUTE  CGRSK'A  PUUSE  "SPECf  R'u!^"  AMPL I  f ijDE  SENSED 
IPL*IPH»  fixed  P3INT  L3W-  AND  KI-PREC  SeJNDS  FSR  INTERP9LATI.9N 

FLjPHs  FL9ATING-S’3I\T  .  - - - 

PLRi°HRs  RUDDER  .  COUPLING  C'^EFFICIENTS  F0.^  I NTERPeLATISN  B9UNDS 
PLF>PME»  ELEVAT3?  .  ,  ,  , 

pLW>pH:'<«.  .VflN.S.. .  .  ■  •  .  .  .  , 

PATIS*  INT.EP^StATnN  SCALER 

_ P3^flAp^.g...C3U=LI'-n  CSEPFICIENT  IMER^eLATED  T9  TEST  FREQUENCY 

G9*'^P,G3*''Ei3“''V=  :3Mp‘>NENT  \eiSF  CURRENT  SPECTRa'l  DENS ITY 
B?Mb  RADIAN  SANDa'ICTh 
SSRT(6‘=»'^) 

SCRi3CE,SCW=  CD.^c^K.p^j-r  SHeRT-CIRCUIT  N8ISE  CURRENT  INDUCED  IN- 
ANiTT'JKA 

_ SC»  Y9TAL  S^?RT»CIRCjIT  N^IEE  CURREN*^  XJ'^  AMPS) 

ENS-.  equivalent  NilSE  PIE'lD  ( VeLTS’/M^TER)  . 

Fhzs  FREwUtf^CY  (!v 

ENFD^s  EQUIVALENT  k'^ISE  FIELD  (IN  DB  3EL®W  1  V8LT/METER) 

C9NSTANTS  A^••''  VARIABLES  PARTICULAR  T8  STREAMER  SECTION 

D'A  r  f  “a  n  t  E  NNA^'”b  I  £  f  A  N  a  i^r'  ^8URCE  ("meters  ) . 

I*"AT.  CHACTER  Alp-JANUMERIC  DESCRIOTISN  BF  STREAMER  MATERIAL 
I'A=  material  C80r.--  1»CAN8PY/  2  =  RAD8ME 

KXa  characteristic  DIME\SI8N  er  dielectric  surface  (METERS) 

STRMIb  streamed  discharge  current  (AMPS) 

XIMs  plsatING-p8!\t  material  C8DE 

X<V»  streamer  spectrum  ce^STAhT  ^ . . . . 

A»  STREAMER  SPECTRJM  CONSTANT 
B»3TREAMER  SPECTRUM  ceNSTANT 
ALP-  streamed  SP-CTRum  CONSTANT 
5ET»  STREAMro  SPECTR^’^  CONSTANT . 

AR3*  STREAME®  SPECTRU*^  TERM 

EXL-  STREAM.fs  spectrum  Term  - - 

GLI.T-  .streamer.. spectrum  TERM.  .  '  . 

DESCRIBED  BELex'-.THE  NBTATISN  IS  AS  FSLLewS 

X-DI3IT  IE  FlDATING  NUMBER  IS  CALLED  EBR . . .  . . 

_ N-DI3IT  IP  FIXED  NUMBER  IS  CALLED  F8R 


•■decimal  PeiNT  (REQUIRED  IN  LSCATISNU  WHEN  SH8WN) 

. . . A. alphanumeric  character  IF  ALPHA  WORD  IS  CALLED  FSR 

E»E  (REQUIRED  WHEN  SHBWN )~  . 

S-SPACE 

+■♦  9R  -  AS  APPR9PIATE 

_ (ALL  F9RMATS  ILLUSTRATED  BEL9W  ASSUME  STARTING  IN  C^BLUMN  1, 

AND  SH8ULD  BE  RIGHT-JUSTIFIED)  ^ -  - 

LA/LANT 

-nsaaaaaaaaaaaaaa  ^  ^ 

■  NNNSS, 

NC9UP  (I3/2X}  ~ 


PSTATllO 

PSTATill 

PSTATUa 

PSTAT113 

PSTATilA 

PSTAT115 

PSTAT116 

PSTAT117 

PSTATlia 

PSTATH9 

PSTATiaO 

_PS.TAJiai 

psTAfiaa' 

psTATiaa 

PSTATia4 

PSTATjaS 

PSTATiaS 

PSTAT127 

PSTATiaS 

PSTATiaS 

PSTAT130 

PSTAT131 

PSTATlSa 

PSTAT133 

PSfAfiiA' 

PSTAT135 

PSTAT136 

PSTAT137 

PSTAT138 

PSTAT1J9 

pstatUo" 

PSTAT141 
PSTAT142 
PSTAT143 
PSTATi44 
PSTAT145 
PSTAT146 
PSTAT147 
PSTAT148 
PSTAT149 
PSTAtiBb 
PSTAT151 
PSTATl 
PSTAT153 
PSTATISA 
PSTAT155 
PSTAT156 
PSTAT157 
“PSTATlSS 
PSTAT159 
PSTAT160 
PSTAT161 
PSTATlBa 
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»  I  J 

n  n|n  o  n  n  o  OjO  n  n  n  o  njn  noon  o  n  o'n  n  o  njn 


_ ESTP<>»iST?>=?ST9  (t9«?ilX»>:S«g<1X>E:9«gjgX) _ 

*+X*XXE*N-'!S+v;xXE+\!^S+X.XXE+N\'SS -  - 

NPJV  (15/2X) 

»K’NKSS  . .  . .  “  . . .  . . 

(II, py) 

■  -NSS .  ■  ■  . . . - . 

_ IT  (6AE<2v;) 

bAAAAAAAA  A4  Atsc~~  - - 

XN/S='?^ALT  lX»f6.1,  iy,rA.l#?)() 

=xx.yxsxyxx«>.5xx»xss . 

'  sNSS 

_ (3(P'5«g#lX),lX),  9R«..  DAFT#WiFUSOI 

'XX.vysyy.yyf  yy.xXES  ■ 

(1 3#  EX) 

3*'-\VES 

PREOIJ  (ES.E/Ex) 

=-!'X.XXE*\i‘''S£ 

_ A  AVT  ^  ?ND>i  (?(-E«E,PX))  ♦  eR...  DIERAT 

«XX  •  vy  3sy  V  .'v  vs;^  ~  - 

ICL^/IC  {I1»1X>'’A2)  <  ly,  IMAT 

s’JSAA.AAAAAAAAAAAA 


^<lGC  ),^.'(lCC)/»(l.:C)iIT(6)/IL(l)ArREQU(90)ALANT{7)A  IC(7 
I  ^'FXcs  1 1  X  E£T9  ( IC  O  *  S'  ST  ?  (1 CG  )  #  RST?  ( 1 OC )  /  I  MAT  (  7 ) 

ImTf^PMATS**  '  - - - - 

39  F95i^'AT(Ay,Ff  .P,^x,4(l?rlo.3,7y) ) 

79  F«RyAt(4A?) 

SO  F9Q''AT(  !3,2y ) 

81  P?R'^AT(F9.2/lx,E9.2,iX,£a‘.g,?y) 

82  '^!)(;‘1AT(!l,2y) 

“51  FfRMAT(6A?/?v5  - - - — - - - — 

84  P';PMAT(rE^p,lX,r6,l/lX,FA,l,py) 

85  F9RMATn(F5'.2vrx)/ixy~ . . .  . 

86  F?RMAT(2(F5.2A2vn 

88  ■■  F9RMAf(E9''*2i2xT  . . . . . . . . . . 

89  FSPMAT(I1,1y,7A2) 

^UD  FSRMAT(lHli?!rX/28KS3I  STATIC  FLECTRICITV . htitZlt/Z/i - 

203  F5RMAT(4(10yA24K»»#»DATA  INPUT  ERR0R**** ) i // ) 

204  F9RMAT(5X/  31HP-STAT1C  M8DEL  EVACUATED  F 9 R~;6 Al79H'7T^RAFTj . . 

.|^g,  F9RMAT(5X^1CHSCALE  SIZE/9X/5HSPEEDx8X/8HALTITUDE/8X*10HCLSUD  TYPE 

206  F8RMAT ( 24Xi 5H ( MPm ) / 9X/ 5H ( KFT ) , / )  ^ - ^ — 

!°!  10X/F4.1/10X#7A2i///) _ 

2ob  FBRMAUbX/llHSTART  FREQ«4  4X#  ICHST&P  FREQ  .#!)X/ 7h£)EuTA-F  > - 

209  F9RMAT{7X/5H(MHZ) j 12X/2(5H(MH7) ,AX W/) 

210  FeRMAT(6X>F6.2/10X#P6.2#8x#F5.2i///) - ^ ^ - T" - - : - 

_211  F&RMAT(5X»8H9ECEIVER/ 10X47HANTENNA</<5Xa5HNSISEa 13X«9HINDUCTIBNj/ 

A 5X#9H9AN0WI0TH/ lOXi 4HAREA^ /#6Xi 5H (XHZ > * 13Xj 6H(M*#2) j/) 


■PSTAT163 
PSTAT164 
PSTATi65 
PSTAT166  ■■■ 
PSTAT161 
>StATr68  ■ 
PSTAT169 
PSTAT170 
PSTAT171 
PStAfiTa 
PSTAT173 
PStAtl74' 
PSTAT175 

■  pstatIT? 

PSTAT177 
PSTAtl78' 
PSTAT179 
PSTAtlSO 
PSTAT181 
‘PSTATfSr' 
PSTAT183 
PSTAT184 
PSTAT185 
PSTAT186 
PSTAT187 
"  PSTATiSl' "■ 
PSTAT189 
’)PSTAT190 
PSTAT191 
PSTAT192 
PSTAT193 

'  pstatT5'5 
PSTAT195 
PSTAT196 
PSTAT197 
PSTAtlSa 
PSTAT199 
■'PSyAT200 
PSTAT201 
PSTAT202" 
PSTAT203 
“PSTAT204 
PSTAT205 
■"PSTAT206 
PSTAT207 

■  PSTATS^ya 
)PSTAT209 
PSTAT210 
PSTAT211 

TgTATSlg 

PSTAT213 

PSTAT214 

IPSTAT215 

PSTAT216 
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^F9RMAT(6y,r5.?^13yiF5«?i///)  _ ^ _  PSTAT217 

L4  F0RMAT(5X/3'mTHE  calculated  CHARGING  CURRENT  IS# IPEIO.S, lX/4HAMPS,PSTAT2i8 

I  ^“paqMAT(i|ji  )  PSTAT519 — 

,  rh.VL _ _  _  ...  PSTAT220 


- Sl?-F5R!^ATni-'^>55XAg6HSRI  F-STATIC  MfPEL  (C8\‘TD)j/>  _  PSTAT221 

219  PSPMAKSXAlSHTHt  PR3BABILITY  ISa  lX/P6.4i  lX#25HTHArfHE  CHARQTnG  CU^TATFeT 
ARREN'Ti /a  8Xj  PQMWlLL  -E  GREATrp  TMAN^l?El0t3>  lx^AHAMPS>//)  PSTAT223 

P  ft  r»  Ail  A  ▼  /  i  ir  w  — -  -  -  -  — .  '■  '  ....I—  T  f  i  V  I  T 


21®  P5RMAT(c(5X,SHF?r5UENCY)A5Xi 13HSHeRT-CIRCUlTA2(5X#10HEQUIVALENT)T7PSTAT22A 

- . .PI£LD  A  EX ) 7Xl  5H  (MHZ )  i  9X>  AH  ( HZ ) ,  1  iX,  PSTAT225 

B6H  (  Ai^PS.)  A  lOV ,  9h  ( y SLTS/M  )  #  6Xa  7H (  DBV/M  )  , / )  PSTA'r2'26 

...|2^.._FeRMAT(5Xy4P3rfITlX^E.  ANTENNA  LOCATED  AT  THE  jAA2i//)  PSTAT227 

223  F9R‘^AT(6X«F6.?j  lCy>P6.2/5X/HHN5N-UNireRf^#///)  PSTAt228 

_ ?21  FeRMAT(?Xf24HALi.  CI5CMAPGES  PERMITTEDi///)  PSTAT229 

722  F5RMAT(5X#27MPUDDER  DISCHARCF  PROHIBITED^///)  FsTAT23Q" 

723  E5R‘'AT(5X/3D-'WINo  TIPS  CISCH4RSF  ®R5HIBITeD////)  PSTAT231 

724  F?RMAT(5X/34MrLEvAteR  TIPS  DISCHARGE  PR9HIBI TEOi /// )  PSTAT232 

72*  F9Rmat(5X,4SHPUDDeR  AND  WINS  TIPS  DISCHARGES  PR9H IB  I TEDi /// )  PSTAT233 

726  F9R‘<AT(5X/46^RUDDeR  AND  ELEVATOR  TI°S  DISCHARGES  PRSHIB I TEDa /// )  PSTAT234 

— Z-!T  TIPS  discharges  PROMTP  TTrn  . /// <  PSTAT235 

ICOl  F5RMAT(EX#33'-i  p?p  STREAMERING  SCCJRING  9N  THE  a7A2)  .  PStAT236' 

1002  FOR^ATISX^iewAND  THE  AnjTENNA  ,rF.2,32H  METERS  AFT  OF  THE  FRONT  OF  PSTAT237 

1  The  a 7A? )  PSTAT238 

1003  F^B^^AT(SX,4pHA^;C'  A  MIMMgh  CHARACTERISTIC  DIMENSION  OF  aF5.2a26H  MPSTAT239 

lETERS  OF  THE  DIELECTRIC  a7AP)  PSTAT240 

.1004  F?PMAT(5Xa5?H  ANp  A  DIELECTRIC  AREA  TO  A/C  FRONTAL  ARFA  RATTR  PP  aPSTAT241 

A^5  •?///)  ■  . . . PstAT24F 

1006  F?Rv.AT(5Xic7WA\D  A  fuselage  DIAMETER  OF  AF5t2/7H  METERS)  PSTAT243 

1027  F5RMAT(5X/3=.HThE  CALCULATED  STREAMERINQ  CURRENT  IS  a1PE8.2a5h  AMPSPSTAT244 


A  a//) 

♦OEFinr  constants 

"  PI»*.C*ATAN(1.Q) 
lERRfO 

E^SiL  ■  (1.0/' (3e,0*Pi)) 


l.CE-09 


SELECT  CORONA  OP  STREAMER  I NG  PROGRAM  OPT  ION 
l«CORaNA  PROGRAM/  ?»STREAMERING  PROGRAM 

READ  82/ NSEC T 

BRANCH  'to' APPROPRTate  PROGRAM 'SECTToN . . 

_ GO  TO  (100/1000)/NSECT 


»»»»  CORONA  DISCHARGE  SECTION  (PROGRAM  OPTION  1)  »»** 
100  CONTINUE  ^ - 

♦  INPUT*  ^  - - - - 

READ  ANTENNA  LOCATION  "  ^ ^ - 

REAP  89/LA< (LANT( J)/J*l/7) _ 

INPUT. NUMBER  OF  COUPLING  COEFFICIENTS  TO  BE  READ 

READ  80j  NCBUP _ _ _ 

READ  IN  THE  -NCOUP-  COUPLING  COEFFICIENTS 


PSTAT240 
13  Sf^-jPSTjATE^l 
PSTA't24F 
PSTAT243 
5h  AMPSPSTAT244 
PSTAT245 
PSTAT246 

_ _ PSTAT247 

PSTAT248 

PSTAT249 

PSTAT250 

PSTAT251 

■PSTAT252 

_ PSTAT253 

PSTAT254 

PSTAT255 

PStAT256 

PSTAT257 

PStAT25'8 

_ PSTAT259 

PSTAT260 
_  PSTAT261 
PSTAT262 

_ PSTAT263 

PSTAT264 
_ PSTAT26S 

PSTAt^W 

_ PSTAT267 

PSTAT268 

PSTAT269 

PSTATI70 
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READ  81/  (ES'^C'(J)/A=  =  Te(J)iR!;Te(j)>J  =  l/NC9UP) 
TERS  eUT  N5K-USEC  ?r  ARRAYS 


KC?=NC'?UP+1 
C9  1  J«NCejlbC>l 
EST9( J)»OtC 
WSTP(J)=C.C 


_ _  J)sG«C _ 

I  CeN^TIKiUE  ^  ^  - 

C  READ  MJYBER  ‘^E  OBtJSRAM  CYCLES 

READ  ZOf 

C  oe  LDP®  Ce‘-TRe!,5  F:v‘^r-RA^•  CYCLES 
D9  ?PS 

_ C  _READ  DISEMARDE  CUE''.Ch  C'=^DF 

REA-5'?c/  U-br'  ■■  ~  “■  . . . . . 

C  READ  aircraft  Tvsr 

"  read  Z3't  (1T(J),  j*i/f ) 

C  READ  A/C  SCALE  SIZE/  S®LED/  ALTI-^UDE 
READ  E*/  VNiSPCiALT 
...  .C  J^EAD  EsrDUEHCY  SEL ‘ CT  '''•iDr 

C  r^TrrvAi^c,  OF 

C  “PDE  ,VF,  -j  -  SEuECTED  ER^TJEKCVS  (U®  T0  9C) 

READ  S?/  M&''rr 
C  TES"^  E?R  ^‘‘DE  SELE'*'' 

!  E  ( ^19D E  D_1  »i:-- ~ j  1 

C  M9DE  ■■."E3»  ■b7”’R'':Af  '"bTRT/Rs’tR'/CRLTA-F  (  is  i-'-MZ)' . 

80c  READ  55/  ES'^RT/ EOT® » EDEL 
3D  TO  PD3 

C  M®DE  /NE*  0/  RE, AD  S~  ERr-uEscvs  Tf  EE  EVAiuAtED 

801  READ  50/  NFR 

_C  READ  r;  ‘iRR  EREObE'.CY  ®SI\T£  (  I^  YHZ) 

read  8. S/  (ErECLfj)/  J^nfIT)  '  . . ’  ■' 

803  CONTINUE 

C  READ  ANTFKiSA  ISDUCTItrS  AREA  AND  RECEIVER  BANDWIDTH 
read  56/  AANT, 

C  READ  CLBUND  TYff  (UcIRRJS/  ?«STRATe  CUMUL’JS/  AsERSNTAL  SNOW) 
READ  89/  ICL®/  <IC(J)*  J=l/7> 

C  ^  — — - 

_C  *INPUT  DATA  ERROR  Cwrc<* 

--  ...- 

IF(NC9UR-10C)  7rC/7?C/8E 
730  lE(MeDEE)  9/ ID/? 

9  lF(NER-?0)  10/10/25 

~V5  1e(I9RE-7)  i1/11/E5 - — - 

II  IP(  A1,T-5C»C)  B/.E/25 

is  iFCNSDEn'-A/A/i"  .  . .  .  . . 

8  DF»rSTP-FSTRT 

'  ■lF<DE) ‘ES/BS/S  ■  . . . . . 

3  IF(DE-EDEL)  B5/Z5/4 _ 

C  AlLBW  ReSM  TB  EXPAND  ERROR  CHECK  ^  ! - 

B5  IERR«1 

5  CONTINUE  ” 

C _ 

C  *PR1NT  INPUT  data#  ^  ~~ — - 


PSTAT271 

PSTAT272 

PSTAT273 

PSTAT274 

PSTAT275 

PSTAT276 

PSTAT277 

PStAt278 

PSTAT279 

PSTAT280 

PSTAT281 

PSTAT282 

PSTAT283 

PSTATaar  ' 

PSTAT285 

PSTAT286 

PSTAT287 

PSTAT288 

PSTAT289 

PStAT290‘~ 

PSTAT291 

PSTAT292 

PSTAT293 

PSTAT29A 

PSTAT295 

PSTAT29^’ 

PSTAT297 

PSTAT298 

PSTAT299 

PSTAT300 

PSTAT301 

PSfATSO?'" 

PSTAT303 

PSTAT30A 

PSTAT305 

PSTAT306 

PSTATjOT 

PSTAT351 

PSTAT309 

PSTAT316 

PSTAT311 

PS7At3i2 

PSTAT313 

PSTAT314 

PSTAT315_ 

PStAt3i6 

PSTAT317 

PStAT3l8 

PSTAT319 

"PSTAtaaO 

PSTAT321 

PSTAT322 

PSTAT323 

PSTAT324 
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y V  ALT/  _  (,IC (j.) /^  J  =  1.a7) 


sOs/ 80^ ■ 


PJ?INT  200  ^ 

. _ IF(  lERR)  2Ql/?0:/2Ql _ 

201  PRINT  203 

.2QI...P_R,,INI.20iL/  (IT(J?/  J»l/6) _ 

PRINT  221/  (LAM(j)/  Jal74)'  . 

PRINT  205 

PRINT  206  - - - 

. PRINT  .?y7».  vNy  =P(;^..alt/  (iC(j)/_  g  =  l 

PRINT  2ca  . '■ 

...  .  PRINT  20S 

IP  I M?DFr )  ■  804,305/ 80V- . 

805  PRINT  210/  rSTRT^FgXPj fDLL 
GS  TO  806 

8pA  PRIN-f  223/  PPESv(l)/  PREag(NPR) 

806  CONTINUE 
PRINT -ail. 

PRINT  212/  OND>‘-/  AA\T 

_ G9  re  (?ll/^l?/"15/714/71E/71A/717j/ 

711  PRINT  721  — 


_ G9  re  (?11/^1?/"1S/714/71E/71A/717I/  lePF 

711  PRINT  721  — 

G9  T9  718 

712  PRINT  722 
65  T9  718 

713  PRINT  723 .  ■  . 

G5  T9  71.8 

714  PRINT  724 
59  T9  718  ‘ 

715  PRINT  725- 
^  G9  T9  718 

716  PRINT  726 
38  T9  718 

"717  PRINT  727  ■ — ^ - ■; - 

,7.18  .C9NTINUE  _ 

IP  ERRSR/  then  A?e2'T'  RJN/""ELSE  C^’^TINUE 
IP(IERR)  27/26/27 

27  PRINT  203  . . . . .  . .  ■■ 

PRINT  216 

--G9  T8  995  - - : - 

..2.6_.„,.p,R.INT  .  2.17_.  _  .  _ 

CPNPUTE  The  TSTAL  'CHARQi’NG"'CU'RREN^^^^  THE'  AiRCRAFf  . 

_ CL9U«PL9AT( TCLg) _ 

SP6fA«<  (-2»354E»09r*(SPD**3)  )  +’T4»S76E-06r*(SPD**2) 
APD 

CHQC"  6*0757E-04*sPDFA»CLeU4XN  ^ - - 

■  IF(CHQC"1«E»03)  7C0/700/701 

700  PRSB»2.0/(CHGC*l.E+06) 

_ 08  T9  702  _ 

701  PReB»2.0E+06/( (CMGC#1»0E+06)*«3)  " 

702  IF(ALT-20»0}  704/704/705 _ • 

704  PR9B«PR9B4CL9U*A1.T/20»0 
G9  T9  706 

705  PR6B*PR8B*CLeu*20#0/ALT  ~ 

706  CONTINUE 

PRINT  214/  CHGC  ^  - - - 


PSTAT326 

. . . . PSTAT327 _ 

PSTAT328 

.  . .  PSTAT32? 

PStAt33b 
PSTAT331 
PSTAT332 
.  .  .  .  PSTAT333 

PSTAT334 
PSTAT335 
PSTAT336 

_ PSTAT337 

PSTAt338 

PSTAT339 

PSTAT340 

PSTAT341 

PSTAT342 

_  PSTAT343 _ 

PSTAT344 

PSTAT345 

PStAT346 

PSTAT347 

PSTAt348 

_  _  PSTAT349  __ 

PSTAT350  " 

PSTAT351 

PSTAT352 

PSTAT353 

PStAt354 

_  _ _ _ PSTAT355 

PSTATSSF 

PSTAT357 

PSTAt358 

PSTAT3S9 

PStAf360 . 

_ PSTAT361 

PSTAT55? 
PSTAT363 
PSTAt364  ■■  ■ 

_  PSTAT365 

+  ( 6  #  65E-b4 )  #SPSTAf 366'"' ■“ 

_ PSTAT367 

PSTAT368 

_  PSTAT369 

PSTAT370 

_ PSTAT371 

PSTAT372 

PSTAT373 

^  TltAta74 - 

.  PSTAT375 

PSTAT376 

_ PSTAT377 

PSTAT378 
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PRINT 

PRINT 


PRSe/CHCaC 


218 


C 

c 

c 


•BEGIN  calculation* 


SCALE 
_ D9 


ceuPLiNG  C‘»t'ri; 
32  J«l#NCeuP 

E(J)eESTe( J)*AANT 


lENTS  BY  INDUCTION  AREA 


32 


C 

C 


WIJI-WSTSCjj-AAr^T 
R(J)«PST?(  jjiAAAf'"' . 

continue. 

scale  C9UPLING  ‘■COEFFIClEN'tSliV 

located  at  or  near  a  given 
SCAEAC«(1.0/yj)**(2»5) 


SCALE  ‘SIZE 

extremity 


UNLESS  ANtENNA" f S ' 


IP(LA)3110i311G/30G? 

3009  GO  T9' <3111,  31 12/31  l'373l  10")  VU 

3111  D9  3120  J»l,NCeUP 

W(  J)»W(  J)*SCAFAC . . . . . 

R( J)»R(U)*SCArAC 

"3 120“'C9NtlNUE  - - 

GO  T0  3114 

3112  03  3121  J»1,NC9U» 

E(J)«E(  J)*£CA.PAC 

R(j)»R'(  j)*scafa;: .  . 

3121  CONTINUE 

■fio 


T9  3114 
3113  D9  3122  j«l,NCet» 
E(  J)»E(U)*SCA'^AC 
W( J)»W( J)*SrAPAC 
3122'  CeMTlNUE . . . . 


GO  Te  3114 

3110  D9‘?l2'3  J»l,NCPUo  — - - - - ^ 

E{J_)»E(J)*Sr.AFAC 

W(j)»W(J)*lCAFAC'  ■  ■'  '  . . .  . 

R(J1»P{J)*SCAFAE 

3123C9NTINUE . .  . .  — . . . . 

3114  CSNTINUE _ _ 

c  SCALE  COMPtJNEN'^  DISCHARGE  CURRENTS 
RU0I"O»182*CW3C 
ELEI»0.364#cmGC 

WINI»0»454#cHGC _ 

C  CALCULATE  C6MP9NEKT  S'^ECTRUM  NaRMALIZERS 
D2R«1«037E-06*SCRT(RUDI  ) 

02E»l*037E-6i*S3RT(EL£l  i - - - 

_D2W»1»037E«06*SCRT(WINI ) 

C  INITIALIZE  FREQUENCY  AND  PRESSURE  PARAMETERS 

_ XCBU»FLBAT(NCeUP)-l«Q _ 

IFIMBDEF)  815#816,815  - 

816  F*FSTRT 

“  qfttttt - - - 

815  LF»1 

F*FREQU(LF)  ^ 

817  CONTINUE  _ .  _ 

EX«EXP(-{  (Alt  ♦  0»002*(ALT**2)  )/2,5t ) ) 


■PSTATSyg 


PSTAT380 

PSTAT381 

PSTAT382 

PSTAT383 

PSTAT384 

PSTAT385 


PSTAT386 


PSTAT387 

PSTAT388 

PSTAT38? 

P'StAf390 


PSTAT391 

pstAt39S 

PSTAT393 

"PSTAT394 


PSTAT395 
PStAt396 
PSTAT397 
PSTAT398 
PSTAT399 
PSTATAbO 
PSTAT401 
PStAf4b2 
PSTAT403 
PSTAT404" 
PSTAT405 
PSTAT406 
PSTATA07 
PSfAt4b8 
PSTAT409 
"  PSTAT410 
PSTAT411 
PSTAT4lg 
PSTAT413 
PSTAt4l4  ■ 
PSTAT415 
~PSTAt416 
PSTAT417 
PSTAT418' 
PSTAT419 
PSTAT420 
PSTAT421 
“PSYAt42S 
PSTAT423 
PSTAT424 
PSTAT425 
PSTAT426 
PSTAT427 
"PSTAt428 
PSTAT.429 
PSTAT430 
PSTAT4ai 
PSTAT432 
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PSTAm4 


A»7«053457E+05*( (760.0*EX)**(-0.25) ) 

XNU»3»83767g'-»03*(  (760«0»EX>  *»(0»»8) ) _ 

BEGIN  FREQUENCY  DEPENDENT  CALCULATISN 
35  C9NTINUE _ _ 

test«xn*f' 

IF (TEST  ■  Xr9jr  36i  36#  38 


38  CALL  9VER 

_ G9  9DS_  _  _  _ . 

36  8NEGA«2tO*PT-»F*i  .bE+06'' .  . . 

.  P.REL«.A#SG^TiYNU/?I  )/SilPt(  (8'^£5A»*2)  4.  (ALPHA**2)) 

DSMP»D2P*=REL 

D9Mr«3gr»pqgt_ _ 

DSMi«»D2N*PREL 

calculate  scaled  CfJPLlM  CSEFFICIENTS 
IFL'IP'iXCTEST) 

.  ...  IFH-irL  *  1  .  .  .  . 

FfFL^A-'dFL) 

FH«FL  4-  1»0  _ _ 

»LR"'’nFL+l) 

?HR»R(  ' 

PLE»E(FL+i) . 

BMraE(IPM+l) 

PLW»W(1FL+1> 

PHV/«W(  IFM-i-l  .) _ _ _ 

RAf!A»(TEST-PL)/(FH-FL) 

PRfPL!^  (PHR-P1.R)*'^ATI.9 
people  4  (PhE"'PuE)*PATI8 
PNtiOL'^  *  (Pt-'W-P1,=.‘)*RATI? 

COMPUTE  Re?T(G(9MESA)) 

3?MR«PR*D3f'R 

G9ME»PE*D‘5ME 

.G9MW«PW*D?MW 

CSNO'jTE  SH9Pt-CTR'C'JIT  N9ISE  CURRENT 
39N»2.04PI.#3N!Dw*in00.0 

S30'-'»SQRT(B“'Y)  . . . 

SCR»G8MP»599M _ 

■  SCE"Q8ME*SB9V! 

_  SCW«G9Ny*S99M_ 

'69  f8r30a7302J3T3#304#30§#3b6i3C7)#I9FF 

302  SCR»SCR/100»0  __  _ 

39  T9  30  A  ■  ■  ■■  "■  •  —  -  —  — . 

303  SCW*SCW/100.0 

Qe  T9  3CS  ^  ^ 

304  SCE»SCE/100»0 _ 

Q6  T8  308 

305  SCR»SCR/100»0 _ _ '  _ 

SCW-SCW/100.0  .  ~  ■ 

68  T8  308  _ 

3b6  SCR-SCR/100.0  ^ 

SCE<8CE/100»0 _ 

68’  T9  308 

307  SCE«SCE/100»0 _ 

SCW«SCW/l00t0 


PSTAT436 
PSTAT437 
PSTAT438 
PSTAT439 
PSTAT440 
PSTAT441 
PSf At442 
PSTAT443 
PStAT444 
PSTAT445 
PSTAT446 
PSTAT447 
PSTAT448 
PSTAT449 
PSTAT450 
PSTAT451 
PSTAT452 
PSTAT4S3 
PSTAT454 
PSTAT455 
PSTAT456 
PSTAT457 
PSfAT458 
PSTAT459 
PSTAT460 
PSTAT461 
PSTAT462 
PSTAT463 


PSTAT464 
PSTAT465 
PSTAT466 
PSTAT467 
PSTAT468 
PSTAT469 
PSTAT4 
PSTAT471 
PStAt472 
PSTAT473 
PSTAt474 
PSTAT475 
PSTAT476 
PSTAT477 
PSTAT478 
PSTAT479 
PSTAT480 
PSTAT481 
PSTAT482 
PSTAT483 
P8TAT484 
P8TAT485 
P8TAT486 
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(59  T9  308 

355T9NtTnU?  - ^ ^ - - — 

C  CeHPUTE  T8_TAL  ^'^6RT-CIRCUJT  KfilSE  CU'^RfNT  _ 

SC«SQ''?T(  (3C'?**2)  +  (SCE**2)  +  (SCW**2)')  ~  . .  . 

__C__C9MPyT_E  equivalent  N'9I5E  PI^UD  _ 

ENF»SC/(eME3Ait'PSlL*A.ANtT  . .  . .  . 

FHZ»E*l»CEf06 

ENFDe»2C .0*AL9G ( EnF ) /2 . 303  ^ - - - 

_C  9UTPUT  PESJLTS 

p^ii;^y-“3>5;  p.rHz/cc^c  . 

C  INICPEMENM  P  ANO  TEST  c’SP  p'EO  ^AKSE  CSHPLETE 

IF'(N9bE~')  'SPC»S2:,32C . 

_  821  F»F+FDEL 

IF(F-PSTP>  35/3-3,40 

820  LP’«1,P+1 
FsF3EaO(L'=’) 

IF(LE-NIFR)  35,33,40 
40  33  t3  999  . . 

C _ _ _ _ _ 

r 

c  ****  StPEA^EPINO  S'CTlOSf  (=P53PA!-  9PTl?\'  2)  **** 
lOOC  CONTINUE 
C  ♦**.  INPUT  *** 

C  . 

■  "c  «?f  X  iT'  r  ?5“'AT.^'5\A'iTrHTnN - - -  - . 

>!£A0  73,  (LANT(  J),  J»l,4) 

C  PPAO  A.IPCRAFT  TVPE 

PEAO  S3,  ( IT( J), j«l,6) 

C  RFAO  A/C  scale  Size,  S^FEC,  AL'^ITJOE 
READ  SA,x\|,  s=d,  ALT 

C  READ  FRE-3U3NCY  SELECT  MSOE  - - 

..  C  M9DE  .EO.d  ■  un1P9R!^  FREQUENCY  INTERVALS  FR9m  FSTRT  TS  FSTP  AT 
C  .  .  intervals  9F  poEL 

C  M9DE  tNEtO  »  J3ER  SClECTEO  PREQUENCIES  (iJP  TB  90 
READ  82,N9DPP 
C  TEST  F?R  M9DE  SELECT 

~  tP(MnOEP)  1802,  1*^01  ^ ^ - 

•E3»0#  read  FSTRT, FSTP, DELTA-F  (IN  MHZ) 

1802  READ  85,FSTpt,fotp',FDpl  . 

09  T9  18C3  .  _ 

'  C  "HBDg'  VNEi'Oi  ■RTa5  >iUYLiE«  ^’E'  ESECUENCl^^^^  lVALUATECi~ . 

1801  READ  80,NPR 

C  f<EAD  IN  NFR  pREDUEN’CY  PCInTs  (IN  MHZ) - - - ^ - 

READ  88,(FPEQu(J),j.1,npr) 

1803  CeNTINUE"  ■  ^  ^ - - - 

C  READ  ANTENNA  INDUCTI3N  AREA  AND  RECEIVER  BANDWITH 

READ  86,AANT,BND>(  '  - - 

C  READ  CLOUD  TYPE  (1«CIRRUS,  I»STRAT9  CUMULUS*  »«FR9NTAL  SNOW) 

READ  B9,ICL8,(li(jJ,J"t,7r  - .  - 

C  READ  IN  CHAR5IN5  MATERIAL  C9DE  AND  MATERIAL 

C Material  caoE  i»wi\dshieloa  2«rao9me  ^ 

READ  89, IM, (!MAT(J),J»1,7)  _ ' 

C  READ  IN  ANTENNA  DISTANCE  (METERS)  AFT  9F  RAOBME  OR  WINDSHIELD 


PSTAT487 


PSTAtASS 

PSTAT489 

PSTATASO 


PSTATA91 

PSfAT492 

PSTAT493 

PSTAT494 

PSTAT495 

PSTAT496 

PSTAT497 


PSTAT498 
PSTAT499 
PSTAT500 
PSTAT501 
PSTAT502 
PSTAT503 
PSTAT504 
PSTAT505 

■■■psfAT5ff5 - 

PSTAT507 

PSTAT508 

PSTAT509 

PSTATSlO 

PSTAT511 

PSTAt5r2“'“ 

BSTAT513 

PSTAT514 

PSTAT515 

PSTAT516 

PSTAT517 

"PSTXTSTF— 

PSTAT519 

PSTAT520 

PSTAT521 

PSTAT522 

PSTAT523 

PSTATSfi 

PSTAT525 

PStAtS26 

PSTAT527 

PSTAT528" 

PSTAT529 

“PSTAT530 - 

PSTAT531 

PSTAT532 

PSTAT533 

~PSTAT534 

PSTAT535 

“PSTAT536 - 

PSTAT537 

~PSTaT838 - 

PSTAT539 

PSTAT540 
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C  and  fuselage  diameter  (METERS) 

_ READ  85>  DA^'Tf  wy.FlJSO! _ _ _ _ 

C  READ  IN  rat  1 9  9^  DIELECTRIC  AREA  TS  AlRCRAFf  p'roNTAL  AREA 
_ READ  «6^0IERAT _ 

c  ■'■■■..■  '  . .  •'  ■“ 

C  **  PRINT  INPUT  data  ** _ 


r 

SRINT  SCO 

PRINT  ?DA/(it(J)i j=li6) 

PRINT  2Elt(_LA\'T(ji,j.l#4) 

oRINT”  iOCl/ (  IMAT(  j)/ jsi/7) 

PRINT  lGG2iCA''T,.  (IMAT(j)j  j=l,7) 


PRINT  ICO?/  -^iX/ (ImAT(J)/  jal/7) 

PRINT  1006/  EyS'-I 
PRINT  lOOA/ClERAT 
PRINT -209 
print  206 

_ PRINT  207/Xv/PP:;>ALT/(IC(J)/Jal/7)  _ . 

print  2C? 

PRINT  209 

IF(MeDEP)  'l?C4/ifc5/lS0A 
ISOE  PRINT  £1C/F3T:;T,FsTP/FDE]. 

G9  TS  i?06 

1804  psinT  223/FPE0U(l)/FREgU(KFR) 


1806  CONTINUE 
PRINT  211 

PRINT  ?1H/ENDN/ AA^•T 
PRINT  217 

C  CSMP'JTE  THE  TSTAL.  CHARGING  CURPENT  T9  THE  AIRCRAFT 
CL9U«PL9AT(ICLe) _ _  ■ 


3PDFA*( ( -P.3S4F-0? )»(SPD**3) )+(4tS76E-06)*(SPD»*2)+6#65E-04*SPD 
CHeC»  6.0757E-04»3P'JFA*CLBU*XN 
IF(CHGC-i»F-03)  1700/1700/1701 

1700  PRS9*2.C/(CHGC*l.E+06) 

G9  TS  17C'2 .  . 

1701  PRe58g.CE*06/(  (ChgC»1«OE4-C6)»*3) 


1702  IF(ALT-2CO)  170u,  1704/1705 

1704  ORfl3*_PR?3*CL'’J*ALT/2C.O 

■  3S'  T9  i7Q6  '  . .  . . . 

1705  PReS«'^99*j:j^y*;2C._0/ALT_ 

1706  eSN'TINUE  ..  . . . 

_ PRINT  2I4/CHGC _ ^ ^ ^ _ _ 

PRINT  219/PR93/CHGC 

C  C9MPUTE  streamer  CHARGING  CURRENT  _  _ 

TEMP«OIERAT*CHGC 

_ 39  T9  (1710/ 1711), IM _ ^ _ ^ _ 

1710  TEMP»TEMP»0.5  ' 

1711  STRMl-TEMP 

PRINT  1027/  STRMI  , 

_ PRINT  218  _ ^ ^ _ . 

C 

C  »*  BEGIN  STREAMER  N9ISE  CALCULATIdN  *» _ 


PSTAT542 

PSTAT543 

PSTAT544 

PSTA_T545 

PSTAf546 

PSTAT547 


PSTAT548 

PSTAT549 

PSTATBSO 

PSTAT551 

PSTAT552 

PSTAT553 


PSIAT554 

PSTAT555 

PSTAT556 

PSTAT557 

PSTAT558 

PSTAT559 

>STAT560 

PSTAT561 

PSTAT562 

PSTAT563 

PSTAT564 

PSTAT565 


PSTAt566 

PSTAT567 

PSTAT568 

PSTAT569 

PSTAT570 

PSTAT571 


PSTAT573 

PSTAt574 

PSTAT575 

PStAt576 

PSTAT577 


PSTAT578 

PSTAT579 

PSTATSSO 

PS.TAT581. 

PStAT582 

PSTAT583 

PSTAT584 

PSTAT585 

PSTAT586 

PSTAT587 

PSTAT588 

PSTAT589 

PSTAT590 

PSTAT591 

PSTAT592 

PSTAT593 

PSTAT594 
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C  .C9NVERT^DIELECTPI^_PARAKETERS  T3  TEET  FR9M  METERS 

: _ FUSDI-FUSDI/0»3076 

C  C9KPUTE  CeUPLIN3  FJKiCTIBN  ^SI 
IF(DAFT)  1712,1713,1712 
1713  PSI»3.0 

sa  ra  i7i7 

l^^l2  39  T9  (1716,1715),  IM  ^ - - 

171_S  psI9NA«1»20E»02/(DAFT«FUSDI  ) 

PSI  *PSiaKA*4  4Kf  '  .  ■  . 

. _ S9_TB  1717  _ _ 

'“1716  P’SiaACA»r(DArT)'*i(-4)  )*o7696+6»6E-b5 

PSI«PSI^NA*AANT 
1717  C9NTINUE  ^ 

C  INITIALIZE  FREGUEK'CY  PARAMPTFRC 

iF(M90EEr  13l'57i3l&ViTl5“' . . .  . . 

1816  F.FStRT 

'  G'a  “la' iEi'7  . . . . . . 

1815  LF»1 

e-«FPE2U(LF)  ^ 

1817  C9NTINUE  _  _  _  _ 

xiMsb.oi .  . .  . .  . . 

Xi<V«l_.27E+C5 
■■  .■■x'MU'«SfRf''i7Tl75E-C9) 

A«0»F97 

fToTSoa  ^  : - 

Al.P»1.67E+0'' 

■9Et»3'.'47E+06"“ .  . . 

C  BEGIN!  FRE^SJEMCY  DEP^^JDENT  CALCULATIQNJ 

1835  a‘'.E3A»2.6*FT*F*l',cE+06 . . . .  . 

C  C9MBUTE  F(X,L) 

ARG»WX*e!"EGA/(2,C*XKV)  ^ 

FXL»B • O^PS ! *PS I » ( 1 •0- ( S IN ( APG ) /ARG ) ) 

C  C9>PUTE'LfffLE”3raYE3Ay . ;  . .  . . 

_ _ TL«J1MEGA**2)*(  (A  +  3)'**?) 

't2»  ( ( Ai3ET+3*ALP)**2)  . 

31»ALP*ALP+0MEGA**2) 

- - 5g^5ET»5ET.,;(&M£.3A>*5) - 

GL I T*  ( T1 +Ta )/(  omega*'#!?) #B  1*32) 

C  C6MBUTE  BIG  G  (9ME5A) 

_ G9M»XNU»XIM*XIn**XKV«XKV»GLIT*FXL/PI _ 

C  C9MPUTE  SH9RT  CIRCOT  CURRENT  (3C) 
39H"2*.0*PI*BNDw*lC00#C 
SBeM.SQRKBa*^) 

RG9M«SQRT(G9M) _ _ _ _ 

SC«SBeM*RG9M 

C  C9MPUTE  EQUIVALENT  NBISE  FIELD _ 

IFIDAFT)  1903,1904,1903 
1903  ENF*SC/(eMEQA*EPSIL*AANT) 

1^04  CeNTINUE 

C  SETUP  8UTPUT  AND  PRINT  RESULTS _ _ 

FHZ«F#1.0E+06 

_ IF  (DAFT)  1900/1901,1900 _ _ 

1901  PRINT  39,  F,FHZi  SC  !  ' 


wm- 

PSTAT597 
PSTAT598 
PSTAT599 
PSTAT600 
PSTAT601 
PSTAT602 
PSTAT603 
PSTAT664 
PSTAT605 
PSTAf606 
PSTAT607 
■PSTAT608 
PSTAT609 
PStAt6ld  ■ 
PSTAT611 
PStAf6l2 
PSTAT613 
PSTAT614 
PSTAT615 
'PSTAt616  ■ 
PSTAT617 
PSTAT6i8 
P$TAT619 
■pSTAT6aO*~ 
PSTAT621 
PSTAT622 
PSTAT623 
PSTAt624  “ 
PSTAT625 
■PSTAT626 
PSTAT627 
PSTAf6'28  ■ 
PSTAT629 
PSTAt636  ■ 
PSTAT631 
TSTAT632 
PSTAT633 
'pSf  AT634  ■■ 
PSTAT635 
PSTAT636 
PSTAT637 

TstAY638 

PSTAT639 

PSTAT640 

PSTAT641 

PSTAT642 

PSTAT643 

PSTAT644 

PSTAT645 

PSTAT646 

PSTAT647 

PSTAT648 


1900  ENFDB«20.0*AL86(ENF)/2.303 

_ _ PRINT  39^F,FHZ^SC.ENF#FNrC9 . . 

C  INCREMENT  F  AND  TEST  ESR  fREGUENCV  RANGE  C9>^PLEtE 

1908  CSNTINUE _ 

IF(M90EF)  1820/1521^1820 
881  F»F*FDEL _ 


IF(F-FSTP)  lS3?ilF35/999 

1880  LF«LF*1  . . . . 

F*FREQU(LF) 

!F(LF-NFR)  1935,1235/999 
999  CONTINUE 

STOP _  - _ ^ _ ■ 

e53d  - — — _ 

_ SUBROUTINE  OVER 

PRINT  1 

1  FORMAT!  45HC3UPLIN3  DATA  NON-EXISTENT  BEYOND  LAST  LISTED//) 
RETURN 


PSTAT650 

-PSIAIill 

PSTAT658 

PSTAT653 

PSTAT654 

PSTAT655 


PSTAT656 

PSTAT657 

PSTATSSS 

PSTAT659 

PSTAT660 

PSTAT661 

PSTAT668 

PSTAT663 

PSTAT66* 

PSTAT665 

PSTAT666 
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